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The Role of Density Functional Theory Methods in the Prediction of Nanostructured 

Gas-Adsorbent Materials 
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Integrated Materials Design Centre, University of New South Wales, Sydney NSW 2052, Australia 

With the advent of new synthesis and large-scale production technologies, nanostructured gas- 
adsorbent materials (GAM) such as carbon nanocomposites and metal-organic frameworks are be¬ 
coming increasingly more influential in our everyday lives. First-principles methods based on density 
functional theory (DFT) have been pivotal in establishing the rational design of GAM, a factor which 
has tremendously boosted their development. However, DFT methods are not perfect and due to 
the stringent accuracy thresholds demanded in modeling of GAM (e.g., exact binding energies to 
within ~ 0.01 eV) these techniques may provide erroneous conclusions in some challenging situations. 
Examples of problematic circumstances include gas-adsorption processes in which both electronic 
long-range exchange and nonlocal correlations are important, and systems where many-body energy 
and Goulomb screening effects cannot be disregarded. In this critical review, we analyze recent 
efforts done in the assessment of the performance of DFT methods in the prediction and under¬ 
standing of GAM. Our inquiry is constrained to the areas of hydrogen storage and carbon capture 
and sequestration, for which we expose a number of unresolved modeling controversies and define 
a set of best practice principles. Also, we identify the subtle problems found in the generalization 
of DFT benchmark conclusions obtained in model cluster systems to real extended materials, and 
discuss effective approaches to circumvent them. The increasing awareness of the strengths and 
imperfections of DFT methods in the simulation of gas-adsorption phenomena should lead in the 
medium term to more precise, and hence even more fruitful, ab initio engineering of GAM. 

PACS numbers: 31.15.-p, 71.15.Mb, 81.05.Zx, 81.05.U-, 81.05.Rm 


I. INTRODUCTION 

A. Rational design of gas-adsorbent materials 

Nanostructured gas-adsorbent materials (GAM) are 
the cornerstones of potentially revolutionary advance¬ 
ments in critical and fast growing technological fields like 
molecular sensing, energy storage and harvesting, and en¬ 
vironmental and sustainability engineering. Their excep¬ 
tional high surface to volume ratio, regular atomic com¬ 
position, tunable reactivity, transport properties, and as¬ 
sembling affinity to form supramolecular systems, have 
permitted the realization of timely and cost-effective ap¬ 
plications like, for instance, the detection and removal 
of toxic substances from water and air, dense storage of 
hydrogen and natural gas in solid state matrices, seques¬ 
tration of carbon dioxide from flue gases generated in 
electricity production plants, design of high-performance 
photovoltaic cells, and enhanced long-lasting operation 
of batteries, to cite just a few examples [IH5]. 

Popular families of nanostructured GAM include 
zeolites, metal oxides nanocrystals (e.g., CaO and 
A-I 2 O 3 ), carbon-based nanomaterials (GN, e.g., nan¬ 
otubes, sheets, met-cars, graphite intercalation com¬ 
pounds, and frameworks of organic pillared graphene), 
metal hydrides nanoparticles (e.g., MgH 2 and LiBH 4 ), 
and covalent- and metal-organic frameworks (GOF and 
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MOF). Of all these species CN and MOF (see Fig. 
stand out as some of the most promising GAM for energy 
and environmmetal applications, particularly to what 
concerns the capture and storage of hydrogen (H 2 ), car¬ 
bon dioxide (CO 2 ), and methane gases [HHTK]. 

Currently reported GAM gas-selectivity and storage 
capacities, however, still remain below the stringent com¬ 
mercial targets set by specialized government bodies and 
agencies. For instance, hydrogen storage systems need 
to achieve an overall capacity of 5.5 wt% hydrogen with 
a volumetric ratio of 40 g/L for competitive vehicle ap¬ 
plications [TB], and capture of the 90 % of the carbon 
dioxide produced in the generation of electricity must be 
reached within less than a 35 % of increase in the fi¬ 
nal costs [niiiH]. The search for optimal gas-adsorption 
processes and GAM, therefore, remains an area of very 
active scientific and technological research. 

A key aspect for potential GAM to be successful is to 
find the optimal chemical compositions and pore topolo¬ 
gies to work under specific thermodynamic conditions. 
The number of possible stoichiometric and structural 
GAM configurations is tremendously large, hence in prac¬ 
tice systematic experimental searches based on trial-error 
strategies turn out to be cumbersome and very inefhcient. 
Rational engineering of gas-adsorbent interactions at the 
atomic scale, represents a key notion to achieve success 
on such a design grand challenge in the short and middle 
term. In this context, computational simulation methods 
emerge as invaluable theoretical tools for the screening 
and rational engineering of auspicious GAM. 
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FIG. 1. (a) Representation of the gas species considered in this review and (b)-(d) some of the most popular families of 
nanostructured gas-adsorbent materials: carbon-based nanomaterials, boron-nitride nanostructures and MOF (see text). 


B. Computational simulation techniques for 
modeling of GAM 

Common simulation techniques in the study of GAM 
can be classified into two major categories: “semi- 
empirical” and “first-principles”. In semi-empirical ap¬ 
proaches, the interactions between atoms are modeled 
with analytical functions, known as force fields or clas¬ 
sical potentials, which are devised to reproduce a cer¬ 
tain amount of experimental data or the results of high- 
accuracy calculations. The inherent simplicity of clas¬ 
sical potentials makes it possible to address the study 
of GAM and gas-adsorption processes considering real¬ 
istic thermodynamic conditions and length/time scales, 
with well-established simulation techniques like, for in¬ 
stance, molecular dynamics and grand canonical Monte 
Garlo uni. With the current computational power and 
algorithm development, key features in GAM which 
are directly comparable to observations (e.g., adsorp¬ 
tion isotherms, elastic properties and diffusion coeffi¬ 
cients [3DH21]) can be computed routinely on a stan¬ 
dard office computer. Also, the relevance of quantum 
nuclear effects in gas-adsorption phenomena can be es¬ 
timated accurately with semi-empirical approaches |25F 
[28] . Nevertheless, in spite of the great versatility of semi- 
empirical methods, classical potentials may present im¬ 
peding transferability issues in certain situations. This 
type of drawbacks is related to the impossibility of mim¬ 


icking the features of the targeted material at conditions 
different from those in which the setup of the correspond¬ 
ing force field was performed [29l [30] . 

In this context, the outputs of first-principles calcu¬ 
lations turn out to be crucial. As the name indicates, 
empirical information is not contained on first-principles 
methods, also known as ab initio. The interactions be¬ 
tween atoms are directly obtained from applying the 
principles of quantum mechanics to the electrons and 
nuclei. Transferability issues, therefore, are totally miss¬ 
ing in ab initio approaches. Examples of first-principles 
techniques include, density functional theory and quan¬ 
tum Monte Carlo, to cite just a few. Although these 
approaches are very accurate, they can be also very de¬ 
manding in terms of computational expense. This cir¬ 
cumstance poses serious difficulties to the study of kinetic 
and thermodynamic effects in extended GAM with ah 
initio techniques. Common acceleration schemes within 
first-principles methods entail the use of pseudopoten¬ 
tials Many materials properties can be predicted 

basing exclusively on the behaviour of the valence elec¬ 
trons, and by employing pseudopotential techniques ex¬ 
plicit treatment of the core electrons is avoided in the 
simulations. Pseudopotentials can actually be the source 
of potential errors, however they make also the simula¬ 
tion of heavy atoms and medium-size systems feasible. 
Fortunately, some strategies can be used to minimize the 
impact of the approximations introduced by pseudopen- 
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tials like, for instance, the projector augmented wave |33j 
and linearized augmented plane wave methods [34]. In 
the present critical review, we will concentrate on an¬ 
alyzing specific aspects of the simulation of GAM with 
first-principles methods, obviating the outcomes of in¬ 
dispensable and physically insightful semi-empirical ap¬ 
proaches. 


C. Capabilities of density functional theory 
methods 

Standard density functional theory (DFT) techniques 
(i.e., based on local and semilocal approximations to 
the electronic exchange-correlation energy, typically LDA 
and GGA) have become the ab initio methods of choice 
in the study and design of nanostructured GAM. These 
techniques have been demonstrated to reproduce with 
notable accuracy the interplay of a wide range of interac¬ 
tions in hundreds-of-atoms systems, while keeping within 
reasonably affordable limits the accompanying computa¬ 
tional expense. However, standard DFT methods present 
some well-known limitations in describing gas-adsorption 
phenomena occurring in low-coordinated atomic envi¬ 
ronments (e.g., the binding of small molecules to sur¬ 
faces and cavities). For instance, due to the local na¬ 
ture of the employed energy functionals standard DFT 
methods cannot reproduce the electronic long-range cor¬ 
relations resulting from instantaneously induced dipole- 
dipole, dipole-quadrupole, quadrupole-quadrupole and 
so on interactions. Another common DFT fault is the 
presence of electronic self-interaction errors, which de¬ 
rive from an imperfect cancellation between the auto¬ 
correlation (i.e., spurious interaction of an electron with 
itself) and exchange energies. This type of errors can al¬ 
ter dramatically the description of charge-transfer com¬ 
plexes, chemical reactions, and electron affinities 
I38j . Fortunately, many developments have been real¬ 
ized during the last two decades which have solved part 
of these drawbacks (e.g., nonlocal and hybrid exchange- 
correlation energy functionals, see next section), permit¬ 
ting so to achieve remarkable agreement between theory 
and experiments. Nevertheless, a number of critical as¬ 
pects remain yet challenging to customary DFT methods 
that could be hindering the rational design of GAM. 


D. DFT challenges in GAM design 

One of such challenges consists in accounting for 
the long-range electronic correlations while simultane¬ 
ously amending the electronic self-interaction errors (see 
Fig. § 1 ). Along with the adsorption of gas molecules 
on surfaces, regions of significant electron depletion and 
accumulation may appear which induce strong electro¬ 
static interactions and the shift of electronic energy lev¬ 
els. In those conditions, one could expect that by adding 
a portion of the exact Hartree-Fock exchange energy to 


the selected DFT functional, in order to lower the elec¬ 
tron self-interaction errors to negligible levels, accurate 
binding energies will follow. Reported evidence, however, 
suggests that accounting for weak dispersion interactions 
in charge-transfer processes may turn out to be also de¬ 
cisive [MlllDI- Unfortunately the causes behind such an 
unexpected result are not yet totally understood, due in 
part to the difficulties encountered in the decomposition 
of DFT energies into fundamental portions. Rational¬ 
ization of this complexity clearly is needed for optimiz¬ 
ing the computational load associated to DFT modeling 
of GAM (i.e., accounting for dispersion interactions nor¬ 
mally requires intensive calculations), and also for ensur¬ 
ing the reliability of prospective and already published 
simulation works on the storage and sequestration of 
gases. 

Another important DFT threat consists in accounting 
adequately for the screening of bare interactions in pe¬ 
riodically extended systems (see Fig. [!>)■ Many-body 
energy and Goulomb screening effects can be equally 
important in physisorption and chemisorption phenom¬ 
ena, and the issues encountered in their description 
are a consequence of the pairwise additivity assumed 
in the construction of most DFT functionals. Essen¬ 
tially, the interaction energy between two atoms remains 
unaltered no matter what medium separates them or 
what collective excitations happen in the material |41F 
H5] . In this context, the adiabatic connection fluctuation- 
dissipation (AGFD) theorem has been exploited to cal¬ 
culate correlation DFT energies that incorporate many- 
body higher-order terms. This is the case of the random 
phase approximation to DFT (RPA-DFT ) and 

the DFT-I-MBD |T71 - [5U| methods (see Sec. jllll T^which at 
present are receiving the highest attention. Nonetheless, 
the development of many-body DFT-based methods are 
still on their infancy and the associated computational 
expenses are elevated (typically ranging from two to four 
orders of magnitude larger than standard DFT). The cor¬ 
responding degree of applicability therefore remains still 
limited. 


E. Inherent complexity of DFT benchmark studies 

Most of ab initio works published to date on the de¬ 
sign of GAM are based on standard DFT calculations (see 
Sec.|^for more details). Standard DFT methods are very 
effective in dealing with large atomic systems and a myr¬ 
iad of user-friendly DFT packages, practically available 
to everyone nowadays, have facilitated their widespread 
use. However, chemical accuracy is generally demanded 
on GAM design (e.g., correct binding energies to within 
~ 0.01 eV) and thus, for the reasons highlighted in the 
previous section, employing standard DFT methods may 
not always be adequate. 

In fact, to carry out DFT benchmark tests on the as¬ 
sessment of GAM is of paramount importance for rigor¬ 
ously establishing acceptable balances between compu- 
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FIG. 2. Schematic representation of the two main challenges 
of customary DFT methods in GAM design: (a) reproduction 
of covalent and dispersion interactions on an equal footing and 
(b) description of Goulomb screening and many-body energy 
interactions in extended systems. 


tational feasibility and predictive reliability. Perform¬ 
ing computational studies of such a type, however, is 
not a trivial task. First, one has to be familiar with 
the technical and foundational aspects of highly accurate 
quantum chemistry methods (QCM) like Mpller-Plesset 
perturbation theory (MP2), the coupled-cluster method 
with single, double and perturbative triple excitations 
[CCSD(T)], and quantum Monte Carlo (QMC), to cite 
just a few. In those approaches the Schrodinger equation 
of the many-electron system of interest is solved directly 
by handling the corresponding wavefunction. Many- 
electron wavefunctions generally are expressed in real 
space as a linear combination of atomic orbitals (LCAO). 
LCAO do not fulfill the conditions of Bloch functions 
and thus the simulation of periodic systems like crys¬ 
tals and surfaces, although possible, it is not straightfor¬ 
ward with them [5Tti5i] . Consequently, most QCM stud¬ 
ies are performed in supercells containing model cluster 
systems where periodic boundary conditions are not ap¬ 
plied. Meanwhile, QCM are inherently complex and the 
associated computational load in general scales poorly 
with the number of electrons (namely as 

Besides the technical intricacies, a conceptual problem 
also arises when performing benchmark studies in model 


cluster systems: on what grounds is correct to generalize 
the so-obtained conclusions to realistic extended materi¬ 
als? To help in answering this question, we show in Fig.[^ 
the partial density of electronic states (pDOS) calculated 
in graphene, i.e., an extended one-atom thick carbon 
surface, decorated with Ca impurities, and a coronene 
molecule doped with a Ca adatom (i.e., Ca@C 24 Hi 2 ) 
using standard DFT methods (in both systems the ra¬ 
tio of Ca impurities to carbon atoms is the same). We 
note that the coronene molecule is generally considered 
as a large enough system in which to carry out compu¬ 
tational accuracy tests. However, as it can be appreci¬ 
ated in the figure the pDOS computed in coronene dif¬ 
fers greatly from the one obtained in the “equivalent” 
graphene-based system (e.g., compare the distribution of 
electronic d states around the Fermi energy level in the 
two cases). Also, charge-transfer calculations based on 
Bader’s theory show that the Ca adatom donates about 
0.08 e~ to the coronene molecule (i.e., weak ionic in¬ 
teraction), whereas in Ca-decorated graphene it supplies 
~ 0.9 e~ to the carbon surface (i.e., strong ionic inter¬ 
action). In fact, quantum confinement effects can intro¬ 
duce important differences in the electronic structure of 
apparently similar extended and cluster systems, thereby 
leading to completely unlike gas-adsorption behaviours in 
the two situations [S3[SB]. In addition to this, long-range 
electronic exchange and nonlocal correlations, which are 
ubiquitous in gas-adsorption phenomena, normally turn 
out to be disguised in model cluster systems. 


An obvious solution to overcome these likely scaling 
problems, affecting most benchmark studies, consists 
in using highly accurate quantum chemistry approaches 
which, alike to DFT, can handle the simulation of peri¬ 
odic systems. All the necessary benchmark and DFT cal¬ 
culations therefore can be undertaken in a same extended 
simulation cell. In this regard, quantum Monte Carlo 
(QMC) emerges as one of the most promising family of 
methods for benchmarking DFT. Actually, the formalism 
of many-electron Bloch functions has been established 
within QMC since long time ago [57] [55] and several im¬ 
plementations of this approach are already available in 
open-source packages like, for instance, CASINO m, 
QWALK [60], and QMCPACK [61]. Also, the balance 
between accuracy and scalability in QMC is among the 
highest of all quantum chemistry methods (e.g., the ac¬ 
companying computational cost climbs as with 

the number of electrons [62]). It is worth mentioning 
that very recent methodological advances have permit¬ 
ted also the implementation of the MP2 and the CCSD 
approaches in a planewave framework |63H65j . The avail¬ 
ability of this last suite of quantum chemistry techniques 
however is rather limited at the moment, although it is 
likely to improve over the next few years |66j . 
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FIG. 3. Partial density of s, p, and d electronic states cal¬ 
culated in a Ca-C 24 Hi 2 molecule {Top), and in Ca-doped 
graphene (Bottom). The ratio of calcium to carbon atoms 
is the same in both systems, namely 1/24. Calculations have 
been performed with standard DFT methods, and the Fermi 
energy level in both systems has been shifted to zero. (This 
figure has been adapted from work. |138| .l 


F. Focus and organization of this review 

The focus of the present critical review is on the assess¬ 
ment of the accuracy of DFT methods in the prediction 
and understanding of nanostructured GAM. Despite that 
the number of DFT studies published so far on GAM ap¬ 
plications is enormous, there is not yet a well-established 
set of best practice principles for meaningful simulation 
of gas-adsorption phenomena and materials. The main 
objective of this review is to assist in dehning this, by 
critically discussing key contributions and unsolved con¬ 
troversies in the field. Also, we identify here the subtle 
problems found in the interpretation of complex bench¬ 
mark studies and propose likely solutions for them. 

Our analysis is constrained to the design of carbon- 
based nanomaterials (GN) and MOF for applications in 
H 2 storage and GO 2 sequestration and capture. The 
main reason justifying this choice is that these types of 
GAM and problems currently are attracting the highest 
attention within the community of computational and en¬ 


vironmental materials scientists [D HT3llb71 - l69] . Moreover, 
by understanding the processes involved in the adsorp¬ 
tion of H 2 and GO 2 molecules on CN and MOF we may 
infer also those affecting many other substances and ma¬ 
terials. For instance, under normal conditions H 2 and 
CO 2 molecules are both linear and consequently cannot 
sustain any permanent dipole moment (e.g., like N 2 ); the 
predominant electrostatic interactions with the solid sor¬ 
bent then are of hydrogen bond and quadrupolar type, 
respectively. Flowever, when CO 2 molecules receive or 
donate electrons they normally get distorted and become 
polar. Thereby, from an electrostatic point of view GO 2 
can turn out to be similar to other important molecules 
like, for instance, water and amonia. On the other hand, 
CN and MOF are structurally and electronically alike to 
boron-nitride nanostructures and covalent organic frame¬ 
works (see Fig. [^, two families of GAM which at present 
are being investigated also very thoroughly [70H74]. For 
the sake of focus, other important mechanisms exploited 
in gas storage and sequestration processes (e.g., gas bind¬ 
ing in bulk chemical form and thermodynamic desta¬ 
bilization of nanocrystals and clusters [ZMZS]) and re¬ 
lated families of materials (e.g., metal oxides and binary 
and ternary hydride compounds [501 - 155] 1. will not be dis¬ 
cussed exhaustively in the present work. 

The organization of this review is as follows. In the 
next section we briefly summarize the binding energy 
targets pursued in the design of H 2 storage and CO 2 
capture and sequestration GAM, together with a short 
description of CN and MOF. Next, we explain the gen¬ 
eralities of DFT methods and provide essential insight 
into its most popular versions. In Sec. IV we review the 
most recent and relevant works done in the assessment 
of the performance of DFT methods in GAM modeling. 
From them, we draw general conclusions on which DFT 
exchange-correlation functionals can be employed safely 
for the simulation of gas-adsorption processes and mate¬ 
rials. Also, we comment on the subtle problems found in 
the generalization of benchmark conclusions obtained in 
model cluster systems to real GAM, and discuss possible 
solutions for them. In Sec. we present a discussion 
on the current tendencies followed in first-principles de¬ 
sign of GAM and propose a number of new directions to 
explore in prospective modeling and benchmark studies. 
Finally, our main conclusions are outlined in Sec. lYg 


II. NANOSTRUCTURED GAM 

Nanostructured gas-adsorbent materials contain char¬ 
acteristic nanoscale motifs which are periodically re¬ 
peated along one, two or three directions. Those motifs 
generally are composed of light atoms which congregate 
to form atomically sparse complexes. There are excellent 
review articles in the literature focusing on the physical 
properties and prospective implementations of GAM (see 
for instance Refs. [551193] 1. hence we highlight here only 
the main traits of some representative species (i.e., CN 
































7 


and MOF). Also, we explain the basic requirements that 
potential GAM must accomplish for achieving effective 
storage and sequestration of H 2 and CO 2 gases. 


A. GAM desiderata 

The “hydrogen storage” and “carbon dioxide capture 
and sequestration” problems have sparked very intense 
research within the communities of chemists, physicists 
and engineers in the last past decade. The discovery of 
new materials with large gas uptake capacities, robust 
thermodynamic stability, fast adsorption-desorption ki¬ 
netics, and affordable production costs, is a key notion 
to succeed in the encountered gas storage and sequestra¬ 
tion challenges 

The binding affinities of potential GAM are deter¬ 
mined by their interactions with the gas molecules, 
which can be of electrostatic, dispersion, and/or orbital 
types [55]. In turn, the strength of the gas-GAM in¬ 
teractions depend on the characteristics of the materials 
and gas species, which in the latter case include elec¬ 
tronic structure, atomic shape, polarizability and perma¬ 
nent dipole and/or quadrupole moments. Key quantities 
in the assessment of gas-adsorbent materials include the 
alignment between frontier molecular orbitals (i.e., the 
energy difference between the highest occupied molecu¬ 
lar orbital -HOMO- in the adsorbate and the lowest un¬ 
occupied molecular orbital -LUMO- in the adsorbent), 
and the resulting charge transfers and binding energies. 
Let us have a closer look into the most relevant physical 
aspects of promising hydrogen and carbon capture GAM. 


1. H 2 storage 

For hydrogen to become the fuel of choice in future 
environmentally clean vehicles and electricity production 
plants, large amounts of H 2 need to be stored within rel¬ 
atively small volumes entailing moderate weights (for a 
detailed list of related targets see Refs. |100[I101) 1. In this 
context, H 2 adsorption on light solid sorbents through 
weak and mild molecule-GAM interactions emerges as 
one of the most promising routes. (Dissociative adsorp¬ 
tion in solid metal hydrides where molecular hydrogen is 
dissociated to form a solid solution with the hydride, is 
also considered as propitious (731751; however, we will 
not discuss such an interesting gas-storage mechanism or 
related family of materials in the present review for the 
sake of focus.) Gurrent estimations of the optimal bind¬ 
ing energy for adsorption of hydrogen at ambient temper¬ 
ature and considering safe delivery pressure conditions 
of 1.5 — 30 bar, amount to ^ 0.2 eV/molecule |102H104j . 
Meanwhile, kinetic effects may be critical for practical ap¬ 
plications (e.g., complete H 2 storage-release cycles must 
be realized within seconds) and when these are taken 
into account favourable binding energies turn out to be 
^ 0.7 eV/molecule [105| . Based on these assessments, the 


binding energy range that typically is targeted in most ab 
initio H 2 -storage works is 0.1 — 1.0 eV/H 2 , which spans 
from dispersion to moderately covalent molecule-GAM 
interactions (see Table I). 

The strength of the H 2 -solid sorbent interactions can 
be tuned by decorating the GAM surfaces with tran¬ 
sition or alkaline-earth metal atoms, and this normally 
tends to increase the affinity towards gas binding |1061 - 
I108j . The orbital mechanism acting behind this effect is 
the Kubas interaction, which consists in electron dona¬ 
tion from the H 2 cr-bonding orbital to the empty metal 
d orbitals and simultaneous electron back-donation from 
the filled metal d orbitals to the H 2 cr-antibonding or¬ 
bital |109l IllOj . A different approach that is promising 
for gas storage at elevated temperatures, is the spillover 
of hydrogen onto inert surfaces. This mechanism im¬ 
plies the chemical activation of H 2 molecules placed on 
top of transition metal sites through catalysis, followed 
by the migration of atomic hydrogen onto the surface of 
the receptor material (e.g., activated carbons, GOF, and 
MOF). Record hydrogen storage capacities of 4 —6 wt. % 
have been recently accomplished with this strategy in 
metal-carbon complexes and MOF [11111115] . 


2. Carbon capture and sequestration 

Carbon capture and sequestration (CCS) actions are 
being implemented in fossil-fuel burning power plants for 
cutting the amounts of green-house gases which are ex¬ 
pelled to the atmosphere. Current CCS means mostly 
rely on the scrubbing of synthesis and flue gases with 
amine solvents. However, solvent-based CCS technolo¬ 
gies generally are not cost-effective due to the high en¬ 
ergy penalty involved in solvent regeneration and equip¬ 
ment corrosion m- Fortunately, membranes and solid 
sorbents do not suffer from these important drawbacks 
and thus they constitute the core of next-generation CCS 
technologies. 

There are two main types of CCS processes in which 
GAM are highly promising: pre-combustion and post¬ 
combustion. In pre-combustion CCS, the fuel is con¬ 
verted into a mixture of hydrogen and carbon diox¬ 
ide gases using processes such as “gasification” or “re¬ 
forming” . The usual thermodynamic conditions in pre¬ 
combustion CCS are high temperatures and high pres¬ 
sures (i.e., above 100° and 20 — 30 bar, respectively), and 
the concentration of CO 2 gas is high. At those condi¬ 
tions, it can be shown that CO 2 sequestration turns out 
to be effective when the strength of the gas-GAM inter¬ 
actions roughly amounts to 2.0 eV/molecule. For this, 
consider the inverse calcination reaction, CaO -I- CO 2 —t 
CaCOa, to be a model of a typical pre-combustion CCS 
process [77H75|. The Gibbs free-energy balance of this 
reaction has been approximated by the function AG = 
— 1.83 -I- 0.0016 • T eV/molecule, where T represents the 
temperature (expressed in Kelvin) |117j . Now, by assum¬ 
ing that the thermal contributions to AG mostly arise 
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FIG. 4. (a) Sketch of a carbon dioxide molecule interacting with a carbon-based surface, (b) Partial density of electronic states 
calculated in the system sketched in (a), expressed in arbitrary units. The electronic states localized in the carbon surface 
are represented in the upper part of the figure, and those corresponding to the gas molecule in the lower, (c) Sketch of a 
carbon dioxide molecule interacting with a Ca-decorated carbon surface, (d) Partial density of electronic states calculated in 
the system sketched in (c), expressed in arbitrary units. The electronic states localized in the Ca-decorated carbon surface are 
represented in the upper part of the figure and those corresponding to the gas molecule in the lower. HOMO bonding and 
LUMO anti-bonding molecular orbitals are indicated in panels (b) and (d) [see text]. 





from the entropy of the gas, ifbind 2.0 eV/molecule 
follows. 

In post-combustion CCS, the CO 2 gas is separated 
from the exhaust of a combustion process. The usual 
thermodynamic conditions in post-combustion CCS are 


low temperatures and low partial CO 2 pressures. In 
that case, it can be shown that the ideal gas-binding 
energy scale for solid sorbents is ~ 0.2 eV/molecule. 
For this, consider that a typical CCS post-combustion 
process can be described by the generic capture reac- 
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tion A + nC 02 —>■ A(C 02 )„, and that the corresponding 
Gibbs free-energy balance vanishes at temperatures close 
to ambient (e.g., 50° degrees Celsius). Now, by assum¬ 
ing that the entropy change in the solid sorbent A upon 
adsorption of molecules is practically null and using tab¬ 
ulated thermodynamic data of the CO 2 gas (i.e., the cor¬ 
responding entropy is ^gas = 5.3 • 10“^ eV/molecule K), 
Abind ~ 0.2 eV/molecule follows. 

Based on these assessments, typical gas-adsorption en¬ 
ergies pursued in most ab initio CCS studies span from 
0.2 to 2.0 eV/C 02 , and depending on the specific aims 
one end or the other of this interval is targeted. As 
it is shown in Table I, that energy range extends from 
weak dispersion to strong covalent molecule-GAM inter¬ 
actions. In addition to these binding energy require¬ 
ments, prospective CCS GAM need to display also fa¬ 
vorable selectivity features with respect to the adsorp¬ 
tion of N 2 , oxygen and water, since those species are also 
abundant in the generated flue and synthesis gases. 

It is worth mentioning that similar electronic orbital 
processes to the Kubas interaction |10911110) explained 
in the previous section, can be exploited also for enhanc¬ 
ing the affinities of GAM towards CO 2 binding. In Fig.|^ 
we represent the partial density of electronic states cal¬ 
culated for a carbon dioxide molecule interacting with a 
pristine carbon surface, and with the same surface deco¬ 
rated with Ga atoms. The presence of calcium dopants 
induces overlappings between electronic p CO 2 orbitals 
and s, d metallic states in the region surrounding the 
Fermi level, resulting in a transfer of charge from the 
molecule to the metal centers. In this process the degen¬ 
eracy of the HOMO liVg bonding and LUMO 27r„ anti¬ 
bonding molecular orbitals is lifted, and electronic charge 
is back-donated from the filled s, d Ca orbitals to the anti¬ 
bonding GO 2 electronic states (compare bottom panels 
in both Figs. and 1^. This orbital mechanism provokes 
the bending of the gas molecule and intensifies the gas- 
GAM interactions |118j . thus providing a route for the 
rational design of sorbent materials for pre-combustion 
GCS applications. 


B. Carbon nanomaterials 

Carbon nanostructures (GN) include an ample range 
of carbon allotropes such as nanotubes, graphene, met- 
cars and graphite intercalation compounds (see Fig. [^. 
GN exhibit a large variety of electronic and transport 
properties resulting from the specific way in which atoms 
are arranged. Carbon nanotubes, for instance, may be 
metallic or semiconducting depending on their diameter 
and the degree of helicity present in their hexagonal-ring 
walls [nil [HO]. Met-cars (i.e., metallocarbohedrynes) 
are extremely stable symmetric clusters formed by metal 
and carbon atoms with stoichiometric formula M„Gm, 
where M typically stands for a transition metal atom 
(M = Ti, V, Zr, Hf, Fe, Cr, and Mo) and n,m = 8,12 
(see Fig. Graphite intercalation compounds 


(not shown in Fig. are composed of alternating planes 
of transition, alkali or alkaline-earth metal and G atoms 
disposed in triangular and hexagonal lattices, respec¬ 
tively. In these materials relative shifts between succes¬ 
sive atomic planes along the out-of-plane direction may 
occur leading to the formation of stacking patterns |123F 

m- 


Interestingly, the electronic and gas-adsorption prop¬ 
erties of GN can be finely tuned by coating their sur¬ 
faces with dopant species, creating defects, or imposing 
mechanical strains |126H131] . Among these, the dop¬ 
ing strategy with metal adatoms has attracted a lot 
of attention because of its technical simplicity and ex¬ 
pected enhancement of the binding affinity towards gas 
mol ecules t hrough the Kubas interaction |109[ II10) (see 
Sec. IIAll |132H135] . However, transition metal atoms 
in carbon surfaces exhibit a strong tendency towards clus¬ 
tering [136H138] and thereby further developments in the 
field of synthetic chemistry addressing this problem are 
awaited. Alternatively, the creation of nanostructured 
networks of defect sites such as vacancies, nanoribbons 
and islands, may increase considerably the number of 
unsaturated carbon sites. This kind of structural modi¬ 
fications generally enhance the interactions with the gas 
molecules and therefore lead to improved gas storage ca¬ 
pacities [13911140) . 


C. Metal-organic frameworks 

Metal-organic frameworks (MOF) are multi¬ 
dimensional nanoporous structures composed of 
metal ions (or clusters) coordinated to rigid organic 
molecules, which are called linkers. The choice of the 
metal ion and linker species completely determines 
the structure and functionality of the resulting MOF. 
Common organic linkers include bidentate carboxylics 
(e.g., HOOC-COOH), tridendate carboxylates (e.g., 
CgHgOe), 1,4-benzenedicarboxyalte (BDC), and azoles 
(e.g., C 2 H 3 N 3 ) molecules. The great freedom with which 
the linkers and metal ions can be chosen and combined 
is reflected in the more than 20, 000 MOF species that 
have been reported in the last two decades (for an 
extensive review on this topic, see Ref. [141) 1. Among 
these, the A-MOF74 {A = Mg and Ni), MOF-5 (e.g., 
Zn 40 tetrahedra linked by BDC organic molecules) and 
A-BTT {X = Ca, Fe, Mn,Cu, Co, Ni, Cr, and Cd, and 
BTT = 1, 2, 5-benzenetristetrazolate) families emerge 
as three of the most investigated compounds. The 
A-MOF74 structure is based on coordinated carboxyl 
and hydroxy groups (i.e., helical A-O-C rods that 
emanate from 6 -coordinated A centers) and its primitive 
unit cell contains 54 atoms. The MOF-5 crystal has a 
cubic symmetry and its conventional and primitive cell 
contain eight and four 0 Zn 40 (BDC )3 formula units, 
respectively (i.e., 424 and 106 atoms). The primary 
building block of X-BTT is a truncated octahedron (i.e., 
a six [X 4 C 1 ]+^ squares and eight BTT ligands structure) 
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which shares its square faces to form a general cubic 
framework, comprising a total of 210 atoms. 

The polarity of the building units and spatial separa¬ 
tion between the organic linkers, affect profoundly the 
binding strength of H 2 and CO 2 molecules and thereby 
constitutes a rationale for the design of MOF-based 
GAM. Also, decorating the MOF surfaces with alkali, 
alkaline-earth and transition metal atoms can enhance 
significantly the adsorption of gas species. Practical re¬ 
alizations of this last type of functionalization, however, 
remain yet technically difficult due to potential atomic 
coalescence |142H144] . Alternatively, modifications of the 
porous frameworks based on the embedment of metal 
nanoparticles have been recently investigated, producing 
very impressive gas-capacity records |145H147] . 

III. DENSITY FUNCTIONAL THEORY 
METHODS 

A. General considerations 

The wave function of a A^-electron system, 
4t(ri, r 2 ,..., Tat), contains all its physical informa¬ 
tion and it is determined by solving the corresponding 
Schrodinger equation. In real materials electrons interact 
through the Couloumb repulsion and are abundant, thus 

is a complex mathematical function that in most of 
the cases is unknown. In 1965, Kohn-Sham developed 
an ingenious theory to effectively calculate the energy 
and related properties of many-electron systems without 
the need of explicitly knowing tk [14811149| . The main 
idea underlying this theory, called density functional 
theory (DFT), is that the exact ground-state energy, E, 
and electron density, n(r), of a many-electron system 
can be determined by solving an effective one-electron 
Schrodinger equation of the form: 

He//^zcr — ^ (1) 

where index i labels different one-electron orbitals and a 
the corresponding spin state. In particular. 

He// = + K,t(r) + I + H.c(r) , (2) 

and 

i<T 

where Vext represents an external field and I 4 c(r) = 
6Exc/Sn{r) is the potential function that results from de¬ 
riving the exchange-correlation energy, Exc, with respect 
to the electron density. 

The exchange-correlation energy has a purely quantum 
mechanical origin and can be defined as the interaction 
energy difference between a many-electron quantum sys¬ 
tem and its classical counterpart. Electrons are indistin¬ 
guishable particles called fermions and the wave function 


describing an ensemble of electrons must change its sign 
when two particles exchange orbital states. This quan¬ 
tum antisymmetry leads to an effective repulsion between 
electrons, called the Pauli repulsion, which helps in low¬ 
ering their total Coulomb energy. Despite E^c represents 
a relatively small fraction of the total energy, this is an 
extremely crucial quantity for all the physical aspects 
of materials and molecules because it acts directly on 
the bonding of atoms. The Exc[n] functional generally 
is unknown and in practice needs to be approximated. 
Actually, this is the only source of fundamental error in 
DFT methods and depending on the approximation that 
is used the resulting approach may turn out to be valid 
or not for describing the systems and phenomena of in¬ 
terest. 

In standard cases Exc[n] is approximated with the ex¬ 
pression 

E-PPro^ln] = I e“r°"(r)n(r)dr , (4) 

where is made to depend on n(r), Vn(r), and/or 

the electronic kinetic energy r(r) = 5 
(see next sections). Actually, the exact form of the 
exchange-correlation energy is readily known and reads 

Exc[n] = J n{r)dr j , (5) 

where na;c(r,r') = n^{r,r') + ni.{r,r') is the exchange- 
correlation hole density at position r' that surrounds an 
electron at position r. Some important constraints on 
na;c(r,r') have already been established. For instance, 
nx{r,E) must be nonpositive everywhere and its space 
integral is equal to —1. Also, the space integral of the 
correlation hole density is zero. These constraints can be 
employed in the construction of approximate, but physi¬ 
cally correct, Exc[n] functionals. 

In the next subsections, we review the most popular 
Exc[n] approximations which are currently employed in 
computational studies of GAM. In Table I, we summa¬ 
rize the main types of bonded and nonbonded atomic 
interactions taking place in condensed matter systems 
and surfaces. In Table II, we outline the adequacy of 
the considered DFT methods in describing those interac¬ 
tions. In Table III, we list some popular computer sim¬ 
ulation packages that can be used to perform standard 
and more advanced DFT calculations. Finally, the rela¬ 
tive degree of accuracy and computational expense asso¬ 
ciated to those DFT approaches are sketched in Fig. 
Additional details on these topics can be found in some 
recent specialized reviews |150H154j . 

B. Local and semi-Local Exc energy functionals 

In local approaches (e.g., local density approxima¬ 
tion -LDA-), Ex;c is approximated with Eq. @ and the 
exchange-correlation energy is taken to the be that of 
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Interaction type —>■ 
Aspect 

; 

Covalent 

Ionic 

Hydrogen bond 

Dispersion 

Many — body 

Spatial decay 

exp [—r^j 

1 

r 

1 

^ iN> 6) 

System dependent 

Energy scale (eV/molecule) 

1.0-0.1 

1.0-0.1 

0.1 

0.1-0.01 

0.1-0.001 


TABLE I. Characteristic traits of usual bonded and nonbonded atomic interactions taking place in condensed matter systems 
and surfaces. 



Random phase approximation 
[ACFDT, empty orbitals] 


Hybrid / Dispersion DFT 
[% of exact E / E "'] 


Meta- or highly-parametrized E 
[t(r)] 


Generalized-gradient approximation 
[Vn(r)] 


Local-density approximation 

[n(r)] 


# Atoms 

~10 


-100 


f 

-1000 


FIG. 5. Representation of common levels of E^c approxima¬ 
tion within density functional theory together with some gen¬ 
eral features. Regions coloured in red indicate “High” and in 
blue “Low”. The typical size of the systems that can be han¬ 
dled with those approaches are indicated on the right margin 
of the figure. 


an uniform electron gas with density n(r), namely ■ 
The exact expression of the functional is known 

numerically from accurate quantum Monte Carlo calcula¬ 
tions |35[ll55j . In order to deal with the nonuniformity of 
realistic many-electron systems, these are normally parti¬ 
tioned into infinitesimal volume elements which are con¬ 
sidered to be locally uniform. In semilocal approaches 
(e.g., generalized gradient approximation -GGA-), E^c is 
approximated also with Eq. (|4|) and is made to 

depend on n(r) and its gradient Vn(r) |15611157j . 

Both local and semilocal approximations satisfy cer¬ 
tain exact Exc constraints (e.g., some exact scaling rela¬ 
tions and the exchange-correlation hole sum rules) and 
can work notably well for systems in which the electronic 
density varies slowly over space (e.g., bulk crystals at 
equilibrium conditions, see Table II). However, this is 
manifestly not the case of GAM which normally contain 
surfaces and pores where u(r) can change very abruptly. 
Moreover, by construction local and semilocal function¬ 
als totally neglect long-range electron correlations, oth¬ 


erwise known as dispersion interactions, which certainly 
are ubiquitous in gas-adsorption phenomena. 


C. Meta-GGA or highly parametrized energy 
functionals 


Meta-GGA functionals are semilocal in nature (i.e., as¬ 
sume the approximate Exc[n] expression in Eq. Q but 
contain an additional degree of elaboration with respect 
to LDA and GGA: the orbital kinetic energy density, 
'''('') = 5 This new ingredient allows 

one to construct functionals which satisfy some addi¬ 
tional exact constraints, like for instance the correct gra¬ 
dient expansion of the exchange energy up to fourth or¬ 
der. Another interesting feature of this family of func¬ 
tionals is that they can be trained to capture the short- 
and middle-range parts of electronic correlation, which in 
some special cases may be enough to describe correctly 
the binding of atoms (see Table II). 

Meta-GGA functionals are versatile and in general 
do not entail significantly larger computational expense 
than standard LDA and GGA. Examples of this fam¬ 
ily include the suite of TPSS functionals due to Perdew, 
Scuseria, and others |15811159j , and the “Minnesota func¬ 
tionals” (MX, with X = 05, 06, 08, 11, and 12) devel¬ 
oped by Truhlar and his group in the University of Min¬ 
nesota m- Meta-GGA functionals in general can pro¬ 
vide accurate lattice constants in solids, surface energies, 
and molecular binding energies, and are already imple¬ 
mented in popular quantum chemistry and DPT packages 
like, for instance, GAUSSIAN m. NWCHEM [T^ . 
SIESTA |163j . and VASP |164j (see Table III). Neverthe¬ 
less, dispersion interactions cannot be reproduced sys¬ 
tematically with these functionals since they lack of ex¬ 
plicit nonlocal correlation kernels (see Sec. Ill E| below). 


D. Hybrid exchange energy functionals 

Hybrid functionals comprise a combination of nonlocal 
exact Hartree-Fock (HE) and local exchange energies, to¬ 
gether with semilocal correlation energies. The propor¬ 
tion in which both nonlocal and local exchange densi- 
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ties are mixed generally relies on empirical rules. The 
popular B3LYP |165) . for instance, takes a 20 % of the 
exact HF exchange energy and the rest from the approx¬ 
imate GGA and LDA functionals. Other well-known hy¬ 
brid functionals are the HSE proposed by Heyd-Scuseria- 
Ernzerhof |166j , PBEO |167j , and the family of Minnesota 
meta hybrid GGA |160| . 

In contrast to local and semilocal functionals, hy¬ 
brids describe to some extent the delocalization of the 
exchange-correlation hole around an electron. This sit¬ 
uation becomes of particular relevance in chemisorption 
and charge transfer processes, where the atomic bonds 
turn out to be elongated or shortened. Also, in strongly 
correlated systems containing d and / electronic orbitals 
life, for instance, transition metal oxides. In other words, 
hybrid functionals are useful in situations where the elec¬ 
tron self-interaction errors, stemming from an imperfect 
cancellation between the artificial interaction of an elec¬ 
tron with itself and the exchange energy, are potentially 
large. Hybrid functionals, however, do not account for 
the long range part of the correlation hole energy and 
thus cannot reproduce dispersion forces (see Table II). Ef¬ 
forts to effectively correct for such drawbacks have been 
made recently by Head-Gordon and other authors (e.g., 
the range separated and dispersion corrected a;B97X-D 
and a;M06-D3 functionals) |168H17l] . 


E. Dispersion-corrected E^c energy functionals 

The condition that any DFT-based dispersion scheme 
must accomplish is to reasonably reproduce the asymp¬ 
totic 1/r® behaviour of the interaction between two par¬ 
ticles separated by a distance r in a gas. Local, semilocal, 
and hybrid energy functionals totally miss this require¬ 
ment. The most straightforward way to correct for such a 
fault consists in adding an energy term to the exchange- 
correlation energy that is attractive and has the form 
Adisp = — (where indexes i and j label differ¬ 

ent particles, and a damping factor is introduced at short 
distances in order to avoid divergences). This approxi¬ 
mation represents the core of a suite of methods termed 
DFT-D which, due to their simplicity and low computa¬ 
tional cost, are being widely used at the moment. Prob¬ 
ably, the most popular family of DFT-D methods are 
the dispersion corrected GGA functionals proposed by 
Grimme m- Despite their appealing features, DFT- 
D methods present several shortcomings. First, many- 
body dispersion effects and faster decaying terms such as 
the Bijjrfj and Cijjrlj interactions are completely dis¬ 
regarded. Second, it is not totally clear from where one 
should obtain the optimal Cij coefficients. And finally, 
once these parameters are determined they remain fixed 
during the simulations, and this can be problematic in 
situations where the orbital hybridization and oxidation 
states change as compared to the free atoms case. 

Several improvements on DFT-D methods have been 
proposed, in which the value of the dispersion coeffi¬ 


cients are made to depend on the specific atomic en¬ 
vironment. Those correspond to the DFT-D3 method 
by Grimme |173j . the vdW(TS) approach by Tkatchenko 
and Schefiler |174] . and the BJ model by Becke and John¬ 
son m- In those approaches the specific variation of 
the Cij parameters are taken on the atomic coordina¬ 
tion or effective volume, and the calculation of reference 
dispersion coefficients and atomic polarizabilities are re¬ 
quired. 

A third degree of elaboration exists in which no ex¬ 
ternal input parameters are needed and the dispersion 
interaction is directly computed from the electron den¬ 
sity. In this context, the exchange-correlation energy is 
expressed as Ej;c = + E^'^^ + A"' where is the 

nonlocal correlation energy. Particularly, A"* is calcu¬ 
lated with a double space integral involving the electron 
density and a two-position integration kernel, (j>, of the 
form 

Ef[n\ = \ JJ n{r')(j){q',q, |r' - r|)n(r)dr'dr . (6) 

In Eq. § , q' and q are the values of an universal func¬ 
tion evaluated in positions r' and r, and ^ is a complex 
function that obeys two main constraints: is strictly 

zero for any system with constant n, and the interaction 
between two molecules has the correct |r' — r|“® depen¬ 
dence for long distances. The described approach was 
introduced by Dion et al. in 2004 m and it repre¬ 
sents a key DFT development since it combines all types 
of interaction ranges within a same formula. Refine¬ 
ments on Dion’s approach have subsequently appeared 
where (i) the original two-position integration kernel (j) 
is modified (e.g., the so-called nonlocal VVIO functional 
due to Vydrov and Voorhis |177Ill79j f. and (ii) the ex¬ 
change term in E^c is substituted with other more accu¬ 
rate functionals (e.g., the so-called vdW-DF2 |180) . vdW- 
optB88 and vdW-optPBE |181j . and vdW-GOOa, |182j 
schemes). All these approaches are termed vdW-DF 
and, thanks to the seminal work of Roman-Perez and 
Soler [183) . their accompanying computational expense 
in a planewave framework is moderate in practice. Nev¬ 
ertheless, the way in which nonlocal correlations are cal¬ 
culated inherently assumes pairwise additivity and thus 
many-body effects are completely disregarded in vdW- 
DF methods (see Table II). 


F. The random phase approximation and 
DFT+MBD 

The adiabatic connection fluctuation-dissipation 
(AGFD) theorem provides a general and exact ex¬ 
pression for the exchange-correlation energy of a 
many-electron system, thereby E^c in principle can 
be calculated in a very accurate way incorporating 
higher-order many-body effects. In particular, the 
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Interaction type — >■ 

DFT flavour 

Covalent 

Ionic 

Hydrogen bond 

Dispersion 

Many — body 

Local and semi — local 
(LDA, GGA) 

13511155H157I 

V 


V/x 

X 

X 

Highly parametrized 
(Meta — GGA, Minnesota) 
|158H160| 



VI X 

VI X 

X 

Hybrid 

(B3LYP, HSE) 

[I65HI63 



V 

X 

X 

Dispersion — corrected 
(DFT - D2, vdW - DF, VVIO) 

rrT^irrTfniTTTirrRnirnrn 


V 

V 

V 

X 

Many — body 

(DFT -t MBD, RPA - DFT) 

|48ll4miT52l[T53l 

V 

V 

V 

V 

V 


TABLE II. Description of the performance of some DFT variants in describing usual types of bonded and nonbonded interactions 
found in condensed matter systems and surfaces. Symbol ^ (x) indicates correct (incomplete) description of the considered 
type of interaction by the corresponding DFT method. 


correlation energy of a system adopts the form 
1 

Ec = - — duj X 

Jo 

X\ (r',r, 


where is the bare density-density response 

function, (r',r,iw) the interacting density-density re¬ 
sponse function at Coulomb coupling strength A, and Tr 
denotes the six-dimensional integration over the variables 
r' and r. The response function x measures the electronic 
response of the system at a point r' due to a frequency- 
dependent electric field perturbation at a point r. In 
the ACFD approach, the adiabatic connection between a 
reference non-interacting system (defined at A = 0) and 
the fully interacting system (defined at A = 1) provides 
the correlation energy of the latter, including many-body 
dispersion as well as other types of electron correlation 
effects. This theoretical framework has been exploited by 
several authors to develop novel many-body DFT-based 
approaches, among which we highlight the random phase 
approximation to DFT and the DFT-I-MBD method. 

In the random phase approximation (RPA-DFT) 
scheme, the interacting response function x\ is defined 
self-consistently via the equation xa = Ao + A |r°-r| ■ This 
approximation has been shown to work reasonably well 
for a number of cluster and extended systems |184| . Fol¬ 
lowing the Adler-Wiser formalism |185[ I186j , xo can be 



iuj) - Xo (r',r,*w) 


(7) 


dX , 


computed by using the occupied and virtual orbitals, and 
the corresponding energies and occupancies obtained in 
DFT calculations. In analogy to post-HF approaches, the 
computational expense associated to RPA-DFT is very 
large (i.e., typically scales with the fourth power of the 
number of particles) and the convergence with respect 
to the basis set generally is too slow |187H191j . Also, it 
must be noted that the short-range part of the electron 
correlation energy is not precisely reproduced by RPA- 
DFT and that this shortcoming may be problematic in 
studying molecular systems m- 

Meanwhile, recent efforts done in the groups of Tkatck- 
enko and Schefller have given birth to the so-called 
DFT-I-MBD method [17Il5()j . which accounts also cor¬ 
rectly for the Coulomb screening and many-body effects 
in electronic systems. In the DFT-I-MBD approach, the 
Schrodinger equation of a set of fluctuating and inter¬ 
acting quantum harmonic oscillators is solved directly 
within the dipole approximation, and the resulting many- 
body energy is coupled to an approximate semilocal DFT 
functional. These approximations result in a signifi¬ 
cant reduction of computational load as compared to the 
RPA-DFT method, allowing one to describe larger sys¬ 
tems containing up to few hundreds of atoms. 

The most appealing features of the RPA-DFT and 
DFT-I-MBD methods is that they are very accurate and 
suitable for studying small-gap and metallic periodic sys¬ 
tems. In fact, these techniques have already been ap- 
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plied successfully to the study of strongly correlated 
crystals and surfaces |193H195) . organic molecules ad¬ 
sorbed on metallic surfaces m, and ionic and semi¬ 
conductor solids |196j . Effective schemes of the RPA- 
DFT and DFT-I-MBD methods have already been im¬ 
plemented in commercial and open-source first-principles 
packages like, for instance, ABINIT |197j . GPAW |203j . 
and VASP |164j (see Table III). 


IV. ASSESSING THE PERFORMANCE OF DFT 
METHODS IN THE DESIGN OF GAM 


Making judicious comparisons between zero- 
temperature calculations and ambient-temperature 
observations turns out to be very complicate due to the 
presence of thermal excitations. Temperature profoundly 
affects the (free) energy and conformation of molecules 
adsorbed on GAM, and the differences with respect 
to T = 0 conditions are difficult to assess with theory 
even when considering the simplest interaction models. 
In fact, free energies and entropies cannot be accessed 
straightforwardly during molecular dynamic simulations, 
in contrast to other thermodynamic quantities like, for 
instance, the total internal energy and pressure. In order 
to evaluate free energies technically and computationally 
involved methods like thermodynamic integration, 
free-energy perturbation and umbrella sampling, have 
to be employed |207H210j . Reasonably then, the most 
direct and exact way of evaluating the performance of 
approximate ground-state methods is to compare them 
with other computational approaches which are known 
to be more accurate. In the case of DFT methods, 
accuracy benchmark tests involve the application of 
quantum chemistry methods (QCM) like Mpller-Plesset 
perturbation theory (MP2), the coupled-cluster method 
with single, double and perturbative triple excitations 
[GGSD(T)], and quantum Monte Garlo (e.g., diffusion 
Monte Garlo -DMC-). 

In what follows, we review recent DFT benchmark 
studies done in the areas of hydrogen storage and carbon 
capture and sequestration problems involving GN and 
MOF. Based on them, we draw general conclusions on 
the suitability of approximate exchange-correlation func¬ 
tionals for the design of GAM. Regrettably, despite that 
the number of these benchmark studies is very modest, as 
it will be shown next, some interpretation and computa¬ 
tional inconsistencies prevail yet among them. Also, we 
discuss the subtle problems found in the generalization 
of DFT performance tests done in model cluster systems 
to extended GAM. A few modeling strategies are pro¬ 
posed to overcome this class of problems, some of which 
are computationally very intensive and some others ten¬ 
tative. 


A. H2 storage 


1. Carbon-based GAM 


A promising alternative to functionalizing carbon- 
based GAM with transition metals atoms is to use al¬ 
kali metal species like lithium and calcium. Aggregation 
issues on alkali metal coatings are expected to be less se¬ 
vere than in the transition metal case (see Sec. IIB), and 
the corresponding H 2 gravimetric densities predicted in 
DFT studies largely surpass the targets set by the U.S. 
Department of Energy (i.e., 5.5 wt. %) |134L I211L1212) . 
Reported experimental records, however, lie significantly 
below the impressive performances anticipated with the¬ 
ory (e.g., a modest 2 wt. %) |213H216) and the reasons for 
those large discrepancies are not yet totally understood. 

In 2009, Cha et al. published a controversial work 
in which the accuracy of standard DFT methods in 
the assessment of H 2 -storage GAM was seriously ques¬ 
tioned m- By using quantum chemistry [i.e., MP2 and 
CCSD(T)] and standard and hybrid DFT methods, Gha 
showed that the binding energies, Fbindj and interatomic 
distances calculated in a model system composed of four 
equidistant hydrogen molecules to a positively charged 
Ca ion [i.e., Ga+( 112 ) 4 , see Fig. [ 61 ], were dramatically 
different. In particular, with standard and hybrid DFT 
functionals favorable molecular binding to the calcium 
cation (i.e., Fbind values around ^ —1.0 and —0.1 eV) 
was found at a H 2 -Ga“'' distance oi z = 2.3 A , whereas 
with MP2 and GGSD(T) methods no effective binding 
was determined at any separation (see Table IV). Gha 
et al. performed analogous gas-adsorption calculations 
in a larger system comprising a coronene molecule (i.e., 
C 24 H 12 ) and a calcium atom, and concluded with the 
same level of inconsistency that found in the Ga+( 112)4 
case. In words of Gha et al. |217] : “(these findings) indi¬ 
cate that previous suggestions for the Ga-based hydrogen 
storage system should be reinvestigated with particular 
care about the charge state of Ga”. The levels of alarm 
associated with this error class, however, were lowered 
shortly afterwards by Ohk et al.. In a formal comment 
on Gha’s work, Ohk argued that the discrepancies be¬ 
tween MP2 and DFT methods reported in the Ga+( 112)4 
and Ga-C 24 lli 2 systems were numerical artifacts stem¬ 
ming from the use of small localized orbitals basis sets, 
i.e., 6 — 311 -b -\-G**, which did not contain polarization 
functions of high enough momenta |218j . By using larger 
basis sets of the Dunning type and performing extrapo¬ 
lation to the complete basis set (GBS) limit to get rid of 
likely finite basis-set errors, Ohk et al. showed that the 
agreement between semilocal DFT and MP2 results ob¬ 
tained in Ga-decorated cluster models was qualitatively 
acceptable (although standard DFT methods exhibited a 
strong tendency towards molecular overbinding, see Ta¬ 
ble IV). The reliability of GGA functionals in the as¬ 
sessment of hydrogen storage materials apparently had it 
been restored (at least, at the qualitative level). However, 
a deep understanding about how standard DFT function- 
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Package 

Periodic/Basis set 

Standard 

Meta - GGA 

Hybrid 

Dispersion — corrected 

Many — body 

ABINIT [I92] 

V/PW 

V 

V 


V 

X 

ADF [198] 

WSTO 

V 


V 

V' (x vdW-DF) 

X 

CASTEP [199] 

^/PW 

V 


y 

y (x vdW-DF) 

X 

GP2K [200] 

V/GT0,PW 

y 


y 

V 

X 

GRYSTAL [M] 

y/GTO 


y 

y 

y (x vdW-DF) 

X 

GAMESS [2ni] 

x/GTO 

y 


y 

y/ (x vdW-DF) 

X 

GAUSSIAN [m] 

Wgto 

V 

y 

V 

y/ (x vdW-DF) 

X 

GPAW [203] 

V/NAO,PW 

V 




X 

MOLPRO [ 20 I] 

x/GTO 

V 

V 

y 

X 

X 

NWGHEM [TH] 

y/GTO,PW 

V 

y 

y 

y/ (x vdW-DF) 

X 

EXPRESSO [205] 

V/PW 

y 


y 

X 

X 

SIESTA [I63] 

^/NAO 

V 


y 

X 

X 

VASP [Ill] 

V/PW 

a/ 



X 

X 

WIEN2K [206] 

y/FP - (L)APW -t lo 




y/ (x vdW-DF) 

X 


TABLE III. Selected computer packages which can be used to perform DFT calculations. “Periodic” refers to the ability of 
handling the simulation of three-dimensional periodic systems. In “Basis set”, “PW” refers to plane waves, “STO” to Slater- 
type orbitals, “GTO” to Gaussian orbitals, “NAO” to numerical atomic orbitals, and “FP-(L)APW-|-lo ” to augmented plane 
waves and local orbitals. Symbol \/ (x) indicates suitability (unsuitableness) of the package to perform a particular calculation 
type. 


als, which definitely incur in self-interaction errors and 
neglect long-range dispersive interactions, could describe 
the adsorption of light molecules on surfaces and cavities 
correctly, was still missing. 

In two subsequent works, the authors of article EH] 
reported further details on the anomalous performance 
of customary DFT methods in describing the fixation 
of H 2 molecules on Ca centers via the Kubas interac¬ 
tion They showed that when a sharp orbital 

transition occurs during adsorption of the gas molecule 
on the metal center, the incompleteness of the electronic 
exchange, present in most DFT functionals, acts criti¬ 
cally by providing erroneous overstabilization of the final 
complex. According to Cha et ai, Ca+( 112)4 exemplifies 
the type of system where the highest occupied molec¬ 
ular orbital (HOMO) changes abruptly from 4s to 3d 
character upon the intake of gas, and the fixation of 
H 2 molecules via the Kubas interaction becomes frus¬ 
trated. Cha’s explanations are based on arguments put 
forward by Gunnarsson more than two decades ago, who 


showed that in situations where the nodal wavefunction 
surfaces are intricate approximate DFT exchange func¬ 
tionals tend to be imprecise by underestimating the en¬ 
ergy cost associated to orbital transitions | 221 j . Cha’s 
arguments indeed brought new physical insight into the 
Ca'*'(H 2)4 contention, however some technical aspects of 
works |219l 1220] could still be criticized (e.g., extrapola¬ 
tion to the CBS limit in the MP2 calculations was not 
pursued). Thus, further benchmark studies were still re¬ 
quired for carefully judging the accuracy of customary 
DFT methods in the design of carbon-based CAM. 

In this regard, Bajdich et al. made an important con¬ 
tribution by performing for the first time quantum Monte 
Carlo (QMC) calculations in the Ca“''(H 2)4 model sys¬ 
tem, and comparing their results to those obtained with 
the MP2, and local, semilocal and hybrid DFT meth¬ 
ods |222) (see Table IV). It was found that QMC calcu¬ 
lations based on the diffusion Monte Carlo (DMC) pre¬ 
dicted no binding at all of the four H 2 molecules to the 
Ca'*' center within the interval 0 < z < 4.6 A , in agree- 
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Work 

Abind (eV/H2) 

-^min ('^) 

dn-H (A) 

Method 

mu 

-1.50 

2.3 

0.81 

LDA/6 - 311 -t -tG * * 


-0.30 

2.3 

0.78 

PBE/6 - 311 -t -fG * * 


-0.40 

2.3 

0.77 

B3LYP/6-311-t-tG** 


-0.05 

4.2 

0.74 

MP2/6-311-t-tG** 


-0.05 

4.2 

0.74 

GCSD(T)/6-311-t-fG** 



-0.90 

-0.20 

-0.15 

2.3 

3.4 

2.3 


PBE/cc - pVQZ 
MP2/cc-pVQZ 
CCSD(T)/GBS 

[2221 

-1.20 

2.3 

0.77 

LDA 


-0.70 

2.3 

0.77 

PBE 


-0.30 

2.3 

0.77 

B3LYP 


No binding 

z < 4.6 

0.77 

MP2 


No binding 

2 < 4.6 

0.77 

HFx - PBEc 


No binding 

2 < 4.6 

0.77 

DMG/ANO - VTZ 


[213] 

-0.16 

3.5 

0.74 

MP2/GBS 


-0.18 

2.2/3.4 

0.77 

AF-QMG/CBS 


TABLE IV. Summary of the binding energy and structural results obtained in the Ca"*"( 112)4 model cluster system by different 
authors employing a variety of DFT and quantum chemistry methods (see text). 
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FIG. 6. Sketch of the (a) Ca^ ( 112)4 model cluster system, 
(b) graphite intercalation compounds | 123H12^ seen from the 
front and top views, (c) organic connector and metal cluster 
model systems, and (d) MOF-5 [also known as IRMOF-1]. 
Cluster systems in (a) and (c) have been widely considered in 
benchmark studies as representative models of GAM depicted 
in (b) and (d), respectively. 


ment with previous MP2 results. In stark contrast to this 
conclusion, popular DFT functionals like LDA, GGA and 
B3LYP vaticinated effective fixation of the gas molecules 
at a distance of 2.3 A . Bajdich’s results are in qualita¬ 
tive agreement (disagreement) with Gha’s (Ohk’s) con¬ 
clusions explained above. Interestingly, the authors of 
work |222] carried out DFT calculations with a blended 
functional consisting of the full HF exchange energy and 
the PBE correlation functional (denoted as HFx-PBEc in 
Table IV). They showed that the results obtained with 
the HFx-PBEc functional were consistent with those ob¬ 
tained with the MP2 and DMC methods (see Table IV). 
In light of those outcomes, Bajdich et al. argued that 
the failure of common DFT functionals in describing the 
Ga'*'(H 2)4 system had its origins on the partial or total 
omission of long-range exchange interactions, in coinci¬ 
dence with Gha’s reasonings found in works |21911220) . 
In spite of the significance of Bajdich’s work, this is nei¬ 
ther free of some technical objections. For instance, in 
the atomic relaxations the H-H intermolecular distances 
were kept fixed to 0.77 A, and the errors stemming from 
the incompleteness of the employed basis sets (i.e., triple- 
zeta) and the fixed-node surface approximation in DMG, 
were not evaluated. 

In a posterior work, Purwanto et al. presented a highly 
accurate study on the binding of the Ga+(H 2)4 complex 
which relied on auxiliary-field QMG calculations 1 ^ . 
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There, extrapolation to the complete basis set limit and a 
better treatment of the sign problem were accomplished. 
Purwanto et al. found that the potential energy curve of 
the four hydrogen molecules exhibits a double-well struc¬ 
ture with almost equal binding minima of ^ —0.18 eV at 
distances 2.2 and 3.4 A. These results are in good agree¬ 
ment with the MP2/CBS calculations performed by the 
same authors and the CCSD(T)/CBS results obtained 
by Ohk et al. [218] (see Table IV). Regarding the dif¬ 
ferences with respect to Bajdich’s work, it was argued 
that these were likely to be originated by the fixed-node 
approximation employed in the DMC calculations. 

Conclusions.- The main conclusion emerging from 
works [21714220112221 [223] is that standard DPT function¬ 
als tend to overestimate the cohesion of the Ca+( 112)4 
complex, and that the likely reason for this flaw is a de¬ 
ficient treatment of the long-range exchange energy. On 
the other hand, highly accurate QCM results obtained 
by different groups on the same system are not fully con¬ 
sistent and therefore a conclusive verdict on the general 
performance of DPT methods, i.e., considering all its pos¬ 
sible flavours, cannot be emitted. The current status of 
computational work on this topic is clearly unsatisfac¬ 
tory. 

Yet, the situation becomes even more puzzling when 
the resemblance between the model system Ca+( 112)4 
and real carbon-based H 2 -storage GAM is brought into 
examination (e.g., see Pig. |^-b). Namely, there is still 
the unresolved question about on which grounds the re¬ 
sults obtained in model cluster systems can be translated 
(if possible at all) to realistic extended materials [55| |S 0 | . 
Por instance, the likely presence of valence s and p elec¬ 
tronic states, coming from the carbon atoms in the sor¬ 
bent, is totally disregarded in the Ca+( 112)4 complex. 
Also, the role of the long-range dispersion interactions, 
which are ubiquitous in gas-adsorption processes, nor¬ 
mally turns out to be underestimated in nano-sized sys¬ 
tems. In this last regard, several analysis on the hydro¬ 
gen storage properties of extended Ca-decorated carbon 
nanomaterials based on dispersion-corrected DPT func¬ 
tionals, have been reported recently. 

Wang et al. have studied the H 2 -storage properties of 
Ca-decorated graphene with the vdW-DP approach due 
to Dion m and compared their results to those ob¬ 
tained with standard local and semilocal methods 
Wang’s calculations show that the H 2 binding ener¬ 
gies obtained with the vdW-DP method lie systemati¬ 
cally below those obtained with GGA-PBE (i.e., stronger 
molecule-GAM interactions by ~ 0.05 eV/molecule), and 
generally above the PDA results (i.e., weaker molecule- 
GAM interactions by ~ 0.05 — 0.10 eV/molecule). (Sim¬ 
ilar conclusions have been attained by other authors in 
alike Ca-decorated carbon nanomaterials [2251122b] . 1 In 
light of these outcomes, it is argued that long-range dis¬ 
persion interactions can be also important in the adsorp¬ 
tion of hydrogen molecules on chemically modified carbon 
nanomaterials, even in the cases where the obtained bind¬ 
ing energies are large (i.e., |i?bind| > 0.1 eV/molecule). 


Actually, Wang’s results appear to expose a new failure of 
standard DPT methods in the assessment of GAM, this 
time related to the description of the electronic correla¬ 
tion energy. Unfortunately, the authors of work [224] did 
not report any comparison with respect to hybrid func¬ 
tionals or QCM calculations, hence the size of the bias in¬ 
curred by the vdW-DP method itself, i.e., due to the ap¬ 
proximations performed on the exchange functional and 
pairwise additivity, cannot be inferred. 

In work EH], Wong et al. have investigated also the 
adsorption of H 2 molecules on metal-decorated graphene 
but considering a large set of transition and alkali metal 
species. The employed method there is an improved ver¬ 
sion of Dion’s approach, the so-called vdW-DP2, in which 
the accuracy of the employed semilocal exchange func¬ 
tional is bettered [180] . In the case of Ca-based coat¬ 
ings, Wong’s results are in qualitative agreement with 
those obtained by Wang et al. [224] however the differ¬ 
ence between the two reported vdW-DP and vdW-DP2 
H 2 -binding energies is not negligible (i.e., about 30 meV). 
When considering other types of dopants, vdW-DP2 
calculations generally predict weaker molecular binding 
than obtained with semilocal GGA methods by ^ 0.1 eV 
per molecule. In fact, Wong’s work comes to reinforce the 
idea that dispersion interactions can affect profoundly 
the interplay between extended materials and hydrogen 
molecules. Gonsequently, nonlocal correlations must be 
taken into account by any DPT functional that is in¬ 
tended for modeling of GAM. 

Conclusions.- Prom all these benchmark studies, we 
can draw the general conclusion that both long-range ex¬ 
change and dispersion electronic interactions are pivotal 
in describing the binding of H 2 molecules to GN-based 
GAM. As for local and semilocal DPT functionals, these 
two elements are totally missing in them, therefore they 
are likely to provide unreliable results on the hydrogen 
storage topic. Situations in which local and semilocal 
GAM predictions seem to be correct normally are for¬ 
tuitous and indicate the presence of large energy error 
cancellations [22511226] . Namely, the sign of the missing 
contributions to the exchange-correlation energy in lo¬ 
cal and semilocal approaches are opposite and therefore 
can compensate each to the other. These error cancella¬ 
tions, however, do not occur systematically and may de¬ 
pend on the specific details of GAM (e.g., see work [227] 
where the sign of the difference between the GGA and 
vdW-DP2 energies varies with the dopant species), and 
thus the use of standard DPT methods for modeling of 
carbon-based H 2 -storage materials is not recommended. 
Meanwhile, new quantum Monte Carlo (QMC) and RPA- 
DPT simulations, this time performed in periodic sys¬ 
tems, are highly desirable for rigorously evaluating the 
performance of meta, hybrid, and dispersion DPT func¬ 
tionals (see works [1521 11951 l228H230j for examples of 
applications of such advanced computational methods to 
simulation of relevant materials). QMC and RPA-DPT 
calculations are also necessary for determining the rel¬ 
evance of many-body energy and Coulomb screening ef- 
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fects on the present class of GAM, which so far have been 
systematically neglected. 


2. MOF 

The interplay between H 2 molecules and MOF (e.g., 
MOF-74,5) are dominated by dispersion interactions and 
the local environment surrounding the binding sites in 
the metal clusters [2dlH2d4) . Consequently, employing 
DFT methods that completely neglect nonlocal correla¬ 
tions (i.e., standard and hybrid flavours) and/or carrying 
out benchmark tests in model cluster systems which are 
too small (see Fig. [^-d) turns out to be inadequate in 
the present case. On the other side, the fact that the 
H 2 -MOF interactions are weak implies that effective gas- 
storage can only be achieved at temperatures well below 
ambient conditions (i.e., T ~ 80 K). Thus, a number of 
strategies have been proposed for increasing the affinity 
of MOF towards hydrogen binding |235j . These include, 
the design and control of porosity EM], functionalization 
of the organic linkers |237j , hydrogen spillover |11411115j , 
introduction of open metal sites in the organic linkers and 
metal clusters |142L I238H241| . and decoration of MOF 
surfaces with alkali metal atoms [242H244] . DFT meth¬ 
ods have been intensively employed in the last two men¬ 
tioned approaches hence we revise next the benchmark 
tests undertaken in those areas. 

It was first argued by Lochan et al. that H 2 molecules 
are attracted by open transition metal and alkali centers 
in MOF through donor-acceptor interactions and electro¬ 
statics ESH]- Based on the outcomes of standard DFT 
calculations, those authors concluded that the strength 
of the H 2 -MOF interactions was within the range of de¬ 
sirable binding for ambient gas storage applications (i.e., 
0.3 — 0.8 eV). Lochan et al. also claimed that semilocal 
DFT functionals could be safely employed in the study 
of H 2 -MOF systems because the corresponding leading 
electronic interactions are strong and predominant over 
dispersion ESH]- 

In two recent studies. Sun et al. have evaluated the ac¬ 
curacy of standard, hybrid and meta DFT functionals in 
the prediction of H 2 adsorption on metal-doped organic 
linker systems |24511246j . They have found that when¬ 
ever transition metal, alkaline-earth or alkali metal atoms 
are used, all DFT, MP2 and CCSD(T) methods provide 
quantitatively similar binding energy results, in satisfac¬ 
tory accordance with Lochan’s findings (see Table V). 
In particular, calculations performed with popular DFT 
functionals like PBE, PW91 and M05-2X are in excel¬ 
lent agreement with gold-standard benchmarks obtained 
with the CCSD(T) method (i.e., equal binding energies 
to within ^ 0.01 eV). Meanwhile, hybrid functionals tend 
to underestimate Ebind slightly by 0.05 — 0.10 eV. Anal¬ 
ogous conclusions have been attained also by Dixit et 
al. in a posterior work done in Li-decorated MOF (see 
Table V) fM7\ . 

A further benchmark test confirming the accuracy of 


DFT methods in describing the interplay between H 2 
molecules and chemically functionalized organic linkers, 
has been recently reported by Jiang et al. |248) . In par¬ 
ticular, Jiang et al. have studied the binding of a hy¬ 
drogen molecule to a small C 4 H 3 Li cluster using com¬ 
mon GGA and hybrid DFT functionals, and the highly 
accurate DMC method. As it can be appreciated in Ta¬ 
ble V, notable agreement between all the considered ap¬ 
proaches in Jiang’s calculations is obtained (i.e., within 
~ 0.01 eV). 


Conclusions.- In light of the benchmark results re¬ 
ported in works |245H248] . we may conclude that cus¬ 
tomary DFT methods appear to perform appropriately 
in the simulation of the hydrogen-storage properties of 
chemically functionalized MOF. The physical reason un¬ 
derlying this favorable outcome is that donor-acceptor 
interactions and electrostatics in this family of GAM are 
dominated by short- and medium-range electron-electron 
exchange and correlations. Nevertheless, a note of cau¬ 
tion must be added here. 


In our compilation of benchmark tests we have realized 
a lack of studies analyzing the performance of dispersion- 
corrected DFT schemes in the simulation of extended 
hydrogen-loaded MOF containing open metal sites. In¬ 
deed, dispersion interactions appear to be secondary in 
the present case however, as we will show in Sec. |IVB 2[ 
when H 2 molecules are replaced by CO 2 this type of inter¬ 
actions turns out to be crucial. The polarizability of the 
H 2 molecule certainly is smaller than that of GO 2 (0.79 
and 2.51 A^, respectively |249j 1. however this still must 
have some effect. Actually, it would be very interesting to 
quantify in which proportion dispersion interactions tend 
to lower the DFT energies reported in Table V (which, 
we remind, have been obtained in model cluster systems). 
Also, it would be highly desirable to perform systematic 
studies on the performance of local, semilocal, and nonlo¬ 
cal DFT approaches in periodic simulation of chemically 
functionalized MOF loaded with hydrogen. In this re¬ 
gard, we would like to mention a recent work by Sumida 
et al. in which it has been shown that neither standard 
nor hybrid DFT functionals can reproduce with accuracy 
the measurements done on the binding of H 2 molecules to 
metal-BTT |25()j . Rather, a range-separated hybrid and 
dispersion corrected DFT functional, i.e., the so-called 
a;B97X-D (see Sec. HID), is found to be necessary for 
a correct interpretation of the experimental findings (see 
for instance Table II in work |25()j L Also, Kong et al. 
have recently found that nonlocal interactions are crucial 
to achieve DFT consistency with respect to the H 2 heats 
of adsorption measured in Zn 2 (BDG) 2 (C 6 Hi 2 N 2 ) |251j . 
Further DFT work on the role of the dispersion inter¬ 
actions in hydrogen storage of chemically functionalized 
MOF, is urgently needed in order to avoid likely modeling 
inconsistencies. 




































19 


Work 

System 



Abind (eV/Ha) 





PW91 

PBE 

B3LYP 

M05 - 2X 

MP2 

CCSD(T) DMC 

[245] H 2 - Ca@H4B2C604 


-0.16 



-0.16 

(CBS) 



Ha - Ca@C8H604 

-0.22 

-0.19 

-0.15 

-0.26 

-0.26 

(CBS) 

-0.24 

(CBS) 


Ha - ScOCsHs 

-0.28 

-0.27 

-0.14 

-0.23 

-0.24 

(CBS) 

-0.23 

(CBS) 


H 2 - Ti@C2H4 

-0.39 

-0.37 

-0.26 

-0.38 

-0.42 

(CBS) 

-0.37 

(CBS) 

[2461 

Ha - Li@C8H604 

-0.15 

-0.13 

-0.10 

-0.16 

-0.16 

(CBS) 

-0.16 

(CBS) 

[247] 

H 2 - Li@CioHio04 

-0.42 

(cc - pVTZ) 

-0.18 

-0.16 

(cc - pVTZ) 

(cc 

-0.17 
; - pVTZ) 


piH] 

Ha - Li@C4H3 

-0.11 

-0.12 

-0.12 

(6- 311G[d,p]) 



-0.14 


TABLE V. Results of different benchmark tests carried out in organic linker MOF models decorated with alkaline earth, alkali 
and transition metal atomic species. The types of localized-orbitals basis sets employed in the calculations and the cases in 
which convergence to the CBS limit is achieved, are indicated within parentheses. 


B. CO2 capture and sequestration 

1. Carbon-based GAM 

The adsorption of CO 2 molecules on carbon-based 
GAM can be of physisorption or chemisorption type de¬ 
pending on whether the material surfaces are smooth, 
contain defects or are chemically functionalized. In 
the case of nondefective surfaces, the interactions with 
the gas molecules are dominated by dispersion forces 
thereby the gas is retained on the adsorbent material 
very weakly |252j . 

In 2006, Xu et al. were the hrst in carrying out 
nonstandard DFT calculations on the binding of CO 2 
molecules to pristine graphene |25dj . By using the hy¬ 
brid ONIOM[B3LYP:DFTB-D] method, they found that 
the corresponding gas-adsorption energy was ifbind = 
—0.03 eV. In Xu’s approach, an hybrid DFT evaluation of 
the interactions between CO 2 and a coronene molecule is 
first performed, and a tight-binding dispersion correction 
is subsequently added in order to account for the presence 
of TT-conjugated interactions in the real material. More 
recently, Umadevi et al. have analyzed the same type of 
problem but employing meta-GGA DFT methods (i.e., 
M05-2X), and surprisingly they have found a large ph¬ 
ysisorption energy of ~ —0.10 eV |254j . The reasons for 


the three-fold discrepancy between Xu’s and Umadevi’s 
results remain unclear to us since QGM benchmark calcu¬ 
lations on the strength and nature of the C02-graphene 
or G02-coronene interactions are practically absent to 
date. We only know of a recent work |255j by Lee et al. 
in which the physisorption energy of GO 2 is calculated 
with the MP2 method. In this case, Abind turns out to 
amount ot —0.09 and —0.13 eV, depending on the relative 
orientation between the gas molecule and carbon plane. 
Nevertheless, the basis set of localized orbitals employed 
in Lee’s study is relatively small (i.e., 6-31G**) and the 
model cluster system in which the MP2 calculations are 
performed contains carbon dangling bonds (in opposition 
to real graphene). 

On the other side, the adsorption of gas molecules 
on carbon nanotubes (CNT) has been thoroughly an¬ 
alyzed by Quihonero et al. with dispersion-corrected 
DFT methods (i.e., B97-D/SVP) |256| . By consider¬ 
ing different types of GNT, diameters, and binding sites, 
Quihonero et al. have concluded that C02-adsorption 
is energetically more favourable in the interior than 
in the exterior of nanotubes, in marked disagreement 
with previous reports |257[ 1258) . Also, they have found 
that the strongest GO 2 -GNT interactions occur in the 
(9, 0) and (5, 5) systems where the computed binding en¬ 
ergy amounts to ^ —0.6 eV. Interestingly, based on a 
symmetry-adapted perturbation theory (SAPT) decom- 
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position of the calculated DFT-D interaction energies, 
the authors of work [256] show that (i) dispersion inter¬ 
actions account for the 80 % of the total attractive forces, 
and (ii) electrostatic interactions resulting from the over¬ 
lap between CO 2 and CNT electronic orbitals, although 
secondary, are yet important. We realize, however, a lack 
of high-level ah initio investigations, e.g., employing the 
MP2 and CCSD(T) methods, in the effects of the CNT 
size and curvature on the binding of CO 2 molecules. This 
type of studies would be highly desirable for a consistent 
evaluation of the performance of DFT methods on this 
topic. Similar investigations to these proposed have al¬ 
ready been performed for methane [259] , a molecule that 
is akin to carbon dioxide in terms of the electric dipole 
and quadrupole moments. It is worth noticing that dur¬ 
ing the preparation of this review we became aware of the 
submission of an article by Smith and Patkowski which 
could fill the mentioned knowledge gap in CO 2 -CNT sys¬ 
tems [260] . 

When the surfaces of the carbon nanomaterials con¬ 
tain reactive defect sites like vacancies, holes (i.e., clus¬ 
ters of vacancies) and edges, the interactions with CO 2 
molecules become more intense. For instance, Cabrera- 
Sanfelix has studied the adsorption of carbon dioxide 
on a defected graphene sheet with standard DFT meth¬ 
ods (i.e., DFT-GGA) and has found that molecular ph- 
ysisorption and chemisorption have an energy cost of 
—0.14 and —1.44 eV, respectively [261] . Similar results 
have been obtained by Liu et al. in an equivalent system 
using an analogous computational approach (i.e., — 0.21 
and —1.72 eV, respectively) [262] . However, the authors 
of this last study propose an equilibrium chemisorption 
configuration that is more symmetric than the one ob¬ 
tained by Gabrera-Sanfelix. In a more recent and tech¬ 
nically exhaustive work. Wood et al. have addressed the 
same kind of gas-adsorption problem by employing both 
dispersion corrected DFT and MP2 methods [263] . Es¬ 
sentially, they find that DFT methods incorporating van 
der Waals forces provide adsorption energies which are 
in good agreement with MP2 results, and that GO 2 ph- 
ysisorption in the edges of a zigzag graphene nanoribbon 
occurs with an energy balance of —0.20 eV. 

Conclusions.- In light of the results reported by 
Gabrera-Sanfelix, Liu, and Wood, we can conclude that 
using standard DFT methods for representing the inter¬ 
actions of GO 2 molecules with defected graphene seems 
to be appropriate. It is not clear to us, however, whether 
the suitability of DFT-GGA methods in this case corre¬ 
sponds to a large cancellation between errors or simply to 
a minor role played by the electronic long-range exchange 
and correlations. 

Meanwhile, some functionalization techniques have 
been proposed for increasing the affinity of carbon-based 
GAM towards GO 2 binding, which is desirable for pre¬ 
combustion applications. Among those we highlight 
the decoration of carbon surfaces with nitrogen and 
light metal atoms, for which a number of experimen¬ 
tal and first-principles computational works have been 


performed [11811264H270] . Mo et al. have recently pre¬ 
sented a computational benchmark study of the adsorp¬ 
tion of carbon dioxide on nitrogen-containing hydrocar¬ 
bon molecules m- Specifically, they have carried out 
extensive GGSD(T)/GBS and dispersion-corrected DFT 
calculations in a G 02 / 2 -methylpyridine (the last with 
formula GgHyN) system. Stabilization of this complex 
occurs through an electron donor ( 2 -methylpyridine)- 
electron acceptor (CO 2 ) mechanism and the attractive 
forces between oxygen and hydrogen atoms. The au¬ 
thors of work [271] show that in order to reproduce the 
Abind gold-standard obtained with the GGSD(T)/GBS 
method (i.e., —0.14 eV), both electronic long-range ex¬ 
change and correlations must be taken into account si¬ 
multaneously. According to Mo’s calculations, BLYP-D3 
is among the PBE-D3, BP86-D3, and TPSS-D3 methods 
(see Sec. HIE) the one peforming the best (i.e., providing 
a binding energy of —0.13 eV). 


Conclusions.- The G 02 / 2 -methylpyridine complex is 
therefore a representative example of a system in which, 
despite of its reduced length, both van der Waals in¬ 
teractions and electronic long-range exchange are simul¬ 
taneously important. Also, we note that “low-cost” 
dispersion-corrected DFT schemes like the BLYP-D3 
one m appear to work remarkably well in this sys¬ 
tem. Indeed, it would be extremely interesting to check 
whether the adsorption of CO 2 (and H 2 ) molecules in 
alike cluster systems could be described correctly also 
with this last type of computationally cost-effective ap¬ 
proaches. 


The accuracy of standard and hybrid DFT functionals 
in describing the interactions of GO 2 molecules with Ga- 
decorated graphene has been recently assessed by Gazorla 
et al. In Gazorla’s work, a comparative study be¬ 

tween DFT and MP2 calculations is presented by follow¬ 
ing an original recipe: instead of adopting the customary 
strategy of steadily increasing the size of polyciclic aro¬ 
matic hydrocarbon (PAH) molecules, the concentration 
of Ca dopants in anthracene (i.e., G 14 H 10 , a relatively 
small PAH) is tuned so as to mimic the partial density of 
electronic valence states in Ga-decorated graphene (see 
Fig. ^). The reason for doing that is to constrain the 
size of the cluster system where to carry out the MP2 cal¬ 
culations as much as possible, while still reproducing the 
main electronic orbital mechanisms occurring in the tar¬ 
geted extended material. In this way, the appearance of 
artificial electronic transitions upon gas-loading are pre¬ 
vented in the model cluster system (e.g., see work [ 220 ] 1 . 
Gazorla’s results show that all considered DFT flavors 
predict equilibrium structures which are very similar to 
that obtained with the MP2 method, and energetically 
favourable GO 2 binding (see Figs. [^-c). Nevertheless, 
the differences between hybrid DFT and MP2 energy es¬ 
timations amount to ~ 0.4 eV, and between standard 
DFT and MP2 to 1.0 — 2.0 eV, with DFT methods pro¬ 
viding always the strongest binding. The origins of the 
observed hybrid DFT and MP2 numerical discrepancies 
(see Fig. Izl^) are rationalized in terms of residual self- 
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interaction errors ]Tn\ . The amount of charge trans¬ 
ferred between the gas molecules and Ca-Ci 4 Hio turns 
out to be the same when computed with either hybrid 
DFT or MP2 methods, and long-range dispersion inter¬ 
actions appear to be secondary in the present system 
(see Table 1 in work |272j 'l. Cazorla et al. therefore 
conclude that the strength of the resulting electrostatic 
interactions, which are dominant, must be equal in the 
two compared cases. On the other side, LDA and PBE 
standard functionals overestimate the transfer of charge 
to the CO 2 molecule by 30 — 40% |272j . In light of 
those outcomes, the use of hybrid DFT functionals is rec¬ 
ommended over that of local and semilocal approaches 
for investigation of the gas-uptake properties of AEM- 
decorated carbon GAM. 

Conclusions.- Cazorla’s prescription for choosing a re¬ 
duced cluster system in which to undertake quantum 
chemistry calculations may be useful for optimizing the 
computational expense associated to DFT benchmark 
tests. Also, for justifying the subsequent generalization 
of the attained conclusions to extended materials [272]. 
Nevertheless, only few selected electronic features of the 
targeted system (e.g., partial density of electronic states 
around the Fermi energy level) can be reproduced at 
a time by playing with doping strategies, and proba¬ 
bly only at a qualitative level. In this context, recent 
methodological progress achieved in the resolution of the 
so-called “inverse band structure-problem of finding an 
atomic configuration with given electronic properties” 
problem [2731 1274] (e.g., genetic algorithm searches), 
could be useful. In particular, those techniques could 
be applied to the design of model cluster systems that 
were replicas of the targeted extended GAM in terms of 
the electronic structure. In that case, besides the chemi¬ 
cal stoichiometry, the atomic structure and shape of the 
simulation cell could be varied until finding a suitable 
system in which to perform all the calculations. (Of 
course, electronic long-range exchange and correlation 
corrections should be added somehow afterwards; this 
could be achieved, for instance, by carrying out addi¬ 
tional calculations in both periodic and cluster systems 
with an adequate first-principles method.) We specu¬ 
late that with such an original effective approach it could 
be possible to address important benchmark controver¬ 
sies affecting the performance of DFT methods in GAM 
modeling (e.g., those explained in the present and previ¬ 
ous sections), while avoiding the computational burden 
and bias introduced by the finite size of model cluster 
systems. 
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FIG. 7. (a) Partial density of electronic states calculated in 
Ca-anthracene, 3Ca-anthracene, and Ca-decorated graphene 
using DFT-PBE. (b) Equilibrium structures obtained for the 
adsorption of CO 2 in 3Ca-anthracene using DFT and MP2 
methods, (c) Adsorption energies for CO 2 in 3Ca-anthracene 
obtained with different geometry optimization and total en¬ 
ergy methods. (Figure adapted from work |272p . 


2. MOF 

The interplay between open-metal site MOF and GO 2 
molecules is dominated by dispersion interactions, which 
represent about the 40 — 60 % of the total gas adsorption 
energy P75l[?75] . Besides dispersion, electrostatic and 
orbital interactions are also important in shaping the gas 


affinity of MOF. In particular, puntual charges localized 
in the unsaturated metal sites polarize the GO 2 molecules 
inducing an electric dipole and a bend distorsion in 
them [2751 l277H279j . When low-energy empty d-levels 
are present a forward donation of electrons from CO 2 to 
the metal atoms occurs, which tends to increase the elec¬ 
trostatic contribution to the gas binding [27511278L1280] . 
This last effect is particularly important in MOF-74 and 
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Work 

MOF — type 

DFT flavour 

Ebind 

AH 

AJ^oxp 

[277] 

Mg - MOF74 

PBE 

-0.230 

-0.189 




PBE -h D2 

-0.439 

-0.399 




vdW - DF 

-0.460 

-0.420 




vdW - DF2 

-0.428 

-0.388 




vdW - PBE 

-0.644 

-0.604 




vdW — optB88 

-0.557 

-0.517 




vdW - C09x 

-0.580 

-0.540 


m 





-0.415 



Mg - DOBDC LDA 

-0.542 

-0.531 


PBE 

-0.228 

-0.209 


DFT - D2 

-0.450 

-0.411 


vdW - DF 

-0.528 

-0.490 


vdW - DF2 

-0.500 

-0.479 


vdW — optB86b 

-0.582 

-0.559 


vdW — optB88 

-0.575 

-0.547 


vdW — optPBE 

-0.608 

-0.593 

[284] 

B3LYP -h D 

-0.430 

-0.393 

[283] 

MP2 : B3LYP + D 

-0.480 

-0.443 

-0.458 ±0. 


[255] 

Ni - DOBDC PBE 

-0.124 

-0.091 


DFT - D2 

-0.361 

-0.317 


vdW - DF 

-0.428 

-0.392 


vdW - DF2 

-0.405 

-0.358 


vdW — optB86b 

-0.504 

-0.473 


vdW — optB88 

-0.496 

-0.447 


vdW — optPBE 

-0.516 

-0.483 

[284] 

B3LYP + D 

-0.403 

-0.368 

[283] 

MP2 : B3LYP + D 

-0.455 

-0.420 

-0.410 ± 0.016 


[283] 

[283] 

Co- 

DOBDC 

LDA 

PBE 

DFT - D2 
vdW - DF 
vdW - DF2 
vdW — optB86b 
vdW — optB88 
vdW — optPBE 

-0.408 

-0.100 

-0.322 

-0.407 

-0.375 

-0.455 

-0.452 

-0.477 

-0.379 

-0.086 

-0.301 

-0.386 

-0.337 

-0.419 

-0.414 

-0.453 

-0.370 ±0.020 

[255] 

Cu- 

HKUST 

LDA 

-0.326 

-0.317 





PBE 

-0.092 

-0.097 





DFT 

-D2 

-0.233 

-0.192 





vdW 

-DF 

-0.283 

-0.241 





vdW 

- DF2 

-0.264 

-0.223 





vdW- 

optB86b 

-0.305 

-0.263 





vdW- 

optPBE 

-0.329 

-0.287 


[255] 







-0.246 ±0.085 


TABLE VI. CO 2 -MOF binding energies, Eibind, and heats of adsorption (T = 300 K), AH, calculated with different DFT 
exchange-correlation functionals. vdW-optB86b, vdW-optB88, and vdW-optPBE represent variants of the original vdW-DF 
and vdW-DF2 nonlocal functionals | 1761I180'| due to Klimes and co-workers |181|. Th e vdW-COOa, variant is due to Cooper |182| . 
and the vdW-PBE one is based on the PBE exchange functional [HZ] (see Sec. |IIIl^ . Experimental values, AiL™®', correspond 
to measured heats of adsorption, and all energies are given in units of eV per molecule. 
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BTT species (see Sec. IIC) containing Ti and V open- 
metal sites, where the cr* antibonding state that results 
from the hybridization of oxygen CO 2 lone pairs and 
metal d\ orbitals remains unoccupied [280] . In the case of 
heavily loaded MOF, quadrupolar CO 2 -CO 2 interactions 
turn out to be also very important. 



FIG. 8. Binding energy of a CO 2 molecule on Ca-BTT com¬ 
puted with the PBE-I-D2 approach, expressed as a function of 
the cutoff distance for the pairwise corrections (adapted from 
Ref. |277p . The CO 2 -MOF system in which calculations are 
performed is represented also in the figure. 

Next, we review recent DFT benchmarking works done 
on the screening of MOF for applications in CO 2 cap¬ 
ture and sequestration. Due to the crucial relevance of 
long-range interactions in CO 2 -MOF systems, we dis¬ 
regard here computational studies considering small or¬ 
ganic molecules as model GAM (see, for instance, arti¬ 
cles [28111282] where the C 02 -pyridine and C 02 -benzene 
interactions are investigated) since proper convergence of 
long-range dispersion forces is achieved within distances 
of at least ~ 10 A . This fact is illustrated in Fig. 
which has been adapted from work |277j , where the DFT 
binding energy of a CO 2 molecule on a MOF containing 
Ca open metal sites is represented as a function of the 
cutoff distance that is employed in the calculation of the 
dispersion interactions. 

In Table VI, we enclose the results of gas binding en¬ 
ergy and enthalpy of adsorption, AiJ, calculations per¬ 
formed by Poloni et al. |277j and Rana et al. [283] in 
different MOF, using standard and dispersion corrected 
DFT functionals. For present benchmarking purposes, 
comparisons with respect to quantum chemistry calcu¬ 
lations performed in cluster-size model systems may re¬ 
sult not meaningful because of the convergence reasons 
explained above, and thus the DFT outcomes are com¬ 
pared with experimental heats of adsorption. Enthalpies 
of adsorption at ambient temperature can be estimated 
with the formula AH = Abind + AAzpe + AEte- Here, 
AEzpe and A Ate represent the zero-point energy (i.e., 
Azpe = {wi} being the corresponding vibra¬ 

tional phonon frequencies) and the thermal energy (i.e.. 


Ate = Avib + E^ot + Atransi + kBT for the CO 2 gas phase 
-where ksT accounts for the energy of an ideal gas-, and 
Ate = Avib for the framework with and without the ad¬ 
sorbate) differences, respectively, between the joint and 
disjoint CO 2 -MOF systems. 

Upon comparison of the computed and measured AH 
values shown in Table VI, we can draw the following two 
conclusions: first, standard DFT methods are far from 
reproducing the experimental heat of adsorption trends, 
with LDA (PBE) presenting large underestimation (over¬ 
estimation) bias; and second, among the two families of 
considered dispersion corrected approaches, i.e., semiem- 
pirical Grimme’s and vdW-A (where X ind icates the 
choice of the A)f^^ functional -see Sec. HIE ), the lat¬ 


ter always provides the better agreement with respect to 
experiments (in particular, the vdW-DF and vdW-DF2 
variants show null AH — AH^^p discrepancies within the 
numerical uncertainties in most of the cases). It is worth 
noticing that the value of the AAzpe and AAte differ¬ 
ences generally amount to less than ^ 50 meV, hence the 
major contribution to AH comes from the Abind term. 
Also we note that the structural traits predicted in CO 2 - 
MOF systems with vdW-A methods, not shown here, 
present overall good agreement with the experiments. 

In Table VI, we also include the AH results obtained 
by Valenzano et al. [284] with an hybrid MP2:B3LYP-D 
approach [28511286] . In that scheme, MP2 calculations 
are carried out first in a cluster model system represent¬ 
ing the adsorption site, and a long-range correction is 
added afterwards. The long-range correction is defined 
as the B3LYP-D energy difference between the extended 
and model cluster systems. (The MP2:B3LYP-D method 
can be likewise understood as starting from a B3LYP- 
D calculation for the full periodic structure and adding 
a high-level correction for the adsorption site.) Valen- 
zano’s MP2:B3LYP-D results obtained in the Mg- and 
Ni-DOBDG systems actually are in the same excellent 
agreement with experimental data than found with the 
vdW-DF and vdW-DF2 methods (i.e., null AH - AH<^^p 
discrepancies within the numerical uncertainties). Also, 
they represent an improvement with respect to primary 
B3LYP-D calculations |284j . On one hand, Valenzano’s 
MP2:B3LYP-D work comes to reinforce the accuracy of 
vdW-DF methods in describing the gas uptake proper¬ 
ties of MOF. On the other hand, it demonstrates the 
reliability and efficiency of the intuitive MP2:B3LYP-D 
approach for undertaking benchmark and GAM design 
studies. 


Conclusions. - In light of the compiled results and anal¬ 
ysis, it can be concluded that the use of standard DFT 
functionals must be avoided in the study of the CO 2 ad¬ 
sorption features of MOF. Meanwhile, the suite of vdW- 
DF methods currently represent the best choice for car¬ 
rying out computational studies of such a type, both in 
terms of computational expense and reliability. In this 
regard, we would like to mention that in order to com¬ 
plete our knowledge of the general performance of DFT 
methods in the simulation of CO 2 -MOF systems it will be 
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highly desirable to consider meta-GGA functionals (e.g., 
M05-2X and M06-2X) in prospective studies. Also, it 
would be very interesting to perform advanced RPA-DFT 
and QMG benchmark calculations in those same systems 
for complementing the comparisons reported against ex¬ 
periments, and determining more precisely the role of the 
different energy contributions to the binding of molecules 
(e.g., many-body dispersion effects). Meanwhile, merg¬ 
ing quantum chemistry methods (e.g., MP2) with peri¬ 
odic DFT simulations that incorporate long-range disper¬ 
sion interactions, appears to be an effective strategy for 
calculating accurate binding energies in situations where 
the selected cluster and extended systems are physically 
alike. 



V. DISCUSSION AND PROSPECTIVE WORK 

From all the results and explanations presented in 
Sec. HV] it can be concluded that standard DFT meth¬ 
ods (i.e., LDA and GGA) generally provide results on 
the fixation of H 2 and CO 2 molecules in GN and MOF 
that are correct only at the qualitative level. Actually, 
in the sole analyzed case of hydrogen binding to chem¬ 
ically functionalized MOF coherent agreement between 
DFT-GGA results and highly accurate quantum chem¬ 
istry calc ulation s has been reported by several authors 
(see Sec. IVA2). However, in that particular case full 
consistency between DFT results obtained in cluster size 
and periodic H 2 -MOF systems is still lacking, and fur¬ 
ther investigations are required for determining the pre¬ 
cise role of dispersion interactions. Therefore, our general 
recommendation is to avoid using standard DFT meth¬ 
ods in first-principles modeling of H 2 -storage and CCS 
materials when pursuing accurate binding energies and 
geometries. In stark contrast to this our advice, most 
of the computational studies published to date on gas- 
adsorption topics have relied heavily on the outcomes 
of DFT-LDA and DFT-GGA simulations (i.e., about the 
80 — 60 % of them in the last 14 years). In Figs. [^and[l0| 
we show an estimation of the total number of DFT works 
reported in the areas of H 2 storage and CCS GAM model¬ 
ing, classified according to the employed DFT functional 
and publication date. 

In the case of hydrogen storage (see Fig. [^, we ap¬ 
preciate an important recent decline in the relative num¬ 
ber of modeling studies performed with LDA and GGA 
functionals, as compared to the trend observed during 
the first decade of the present century (i.e., a decrease 
of the ~ 20 %). On the other side, the relative num¬ 
ber of DFT works based on dispersion corrected schemes 
has increased significantly over the same period of time 
(i.e., from ^ 1 % to ^ 10 %). This last datum reflects, 
on one hand, the improved feasibility of electronic band 
structure methods incorporating van der Waals interac¬ 
tions (and the increasing availability through commercial 
packages) and, on the other, the awareness of the impor¬ 
tance of this type of interactions within the community 


FIG. 9. Number of DFT-based “H 2 storage” articles pub¬ 
lished from 2000 onwards, classified according to the employed 
i?a;c approximation (our estimation relies on data extracted 
from the “Web of Knowledge”, December 2014). The solid 
dots represent the total number of works published in different 
years intervals, and the boxes the percentage corresponding 
to each DFT variant. 


of materials scientists. Also, we notice that the popu¬ 
larity of hybrid and meta-GGA functionals have grown 
considerably within the last ten years. Nevertheless, all 
these trends refer to percentages and the truth is that, 
since the total number of published articles has grown 
almost linearly since year 2000, the total number of com¬ 
putational works that are prone to revision (i.e., those 
performed with standard DFT methods) has actually in¬ 
creased during the last years. 


As for the modeling of GCS GAM (see Fig. 10), we 
also acknowledge a steady increase in the total number 
of published works and a significant recession in the per¬ 
centage of recent LDA and GGA studies (i.e., a decrease 
from ~ 60 % to ~ 40 % in the last four years). Mean¬ 
while, the number of dispersion corrected DFT works 
has increased remarkably in the last years, reaching a 
maximum peak of 25 % recently, while the use of hybrid 
functionals has been maintained more or less constant 
around 30 %. Also, meta-GGA functionals are becom¬ 
ing increasingly more popular although these are still the 
least preferred among all the considered DFT variants. 

The ultimate DFT tendencies revealed in both H 2 stor¬ 
age and GGS GAM modeling fields are quite similar, 
namely there is a firm surge in the use of hybrid and 
dispersion corrected approaches in detriment to those of 
LDA and GGA. This fact shows that the outcomes of 
complex DFT benchmark studies are being assimilated 
progressively by the community of materials scientists 
specialized in GAM applications. In spite of this positive 
conclusion, the total number of recent GAM design works 
that still rely exclusively on LDA and GGA calculations 
is, in our opinion, unjustifiably too large. 

In the Introduction section of this review we com¬ 
mented on the two main threats that customary DFT 
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FIG. 10. Number of DFT-based “Carbon capture and stor¬ 
age” articles published from 2000 onwards, classified accord¬ 
ing to the employed E^c approximation (our estimation relies 
on data extracted from the “Web of Knowledge”, December 
2014). The solid dots represent the total number of works 
published in different years intervals, and the boxes the per¬ 
centage corresponding to each DFT variant. 


approaches pose to the design of GAM, namely the dif¬ 
ficulties in (i) accounting simultaneously for the long 
range electron-electron exchange and correlations, both 
of which are omnipresent in gas-adsorption phenom¬ 
ena, and (ii) reproducing many-body energy effects and 
Coulomb screening in extended systems. 


Regarding issue (i), we have demonstrated in Sec. IV 


that good agreement with respect to quantum chemistry 
calculations and experimental data is achieved when well- 
balanced exchange-correlation functionals are employed 
that include van der Waals interactions, on one hand, 
and correct to some extent for the inherent electron self¬ 
interaction errors, on the other. This is the case, for 
instance, of the mixed BLYP-D3 method which in the 
C02/2-methylpyridine system has been shown to per¬ 
form at t he same level of accuracy than CCSD(T)/CBS 
(see Sec. IVB1 and work [271j ). Also, the so-called 
a;B97X-D approach, based on a range-separated hybrid 
and dispersion corrected DFT functional (see Sec. HID), 
has been shown to reproduce correctly the binding en¬ 
ergy trends measured for H 2 molecules on metal-BTT 
complexes (see Sec. IV A 2 and work [250] ). We con¬ 
sequently argue that the safest options for undertaking 
first-principles computational work on gas-adsorption ap¬ 
plications and the design of GAM are to use DFT func¬ 
tionals which incorporate, either in an exact or approxi¬ 
mate way, long range electron-electron exchange and cor¬ 
relations. The performance of these “full long range cor¬ 
rected” (FLRC) DFT functionals has been tested very 
recently by several authors in standard benchmark data 
sets, and very promising results have been attained in all 
the cases (e.g., see studies |I68L 11691117II [T78] ). 

An interesting aspect of FLRG functionals is that these 
can be constructed in principle by taking any pure or long 


range corrected hybrid DFT functional as a start, and 
subsequently adding the missing nonlocality of the cor¬ 
relation energy in the form of additional terms. Also, 
the degree of sophistication and computational expense 
associated to FLRG functionals can be chosen almost 
at wish by adopting simpler or more complex exchange- 
correlation correction schemes. For instance, a computa¬ 
tionally inexpensive FLRG solution may consist in com¬ 
bining the semiempirical dispersion correction approach 
due to Grimme with any hybrid or meta DFT func¬ 
tional of one’s personal taste (e.g., the already employed 
WB97X-D, TPSS-D3, M06-D3, and BLYP-D3 function¬ 
als). A superior FLRG blend, both in terms of accuracy 
and computational expense, could be achieved by con¬ 
sidering full nonlocal functionals like vdW-DF or VVIO 
rather than semiempirical correlation correction schemes 
(see for instance work [128]). In conclusion, we strongly 
recommend to use FLRG functionals in future ab ini¬ 
tio modeling of GAM and to analyze their performance 
comprehensively in prospective DFT benchmark studies. 


Regarding the second customary DFT challenge (ii) 
mentioned above and explained in Sec. |I D[ many- 
body energy and Coulomb screening effects can vary 
considerably the polarizabilities, and consequently the 
forces, of gas molecules interacting with extended sur¬ 
faces and other molecules mm- The RPA-DFT and 
DFT-I-MBD methods emerge as the two DFT variants 
which can deal efficiently with these types of many-body 
effects (see Sec. IIIFl. In particular, the DFT-I-MBD 
method is very well suited for studying large periodic 
systems containing up to few hundreds of atoms [481450] . 
Nevertheless, to the best of our knowledge, neither the 
RPA-DFT nor the DFT-I-MBD method has been applied 
yet to the simulation of hydrogen storage or carbon cap¬ 
ture and sequestration problems, as shown by the lack of 
such studies in the comprehensive analysis presented in 
Sec. IV Prospective RPA-DFT and DFT-I-MBD works 
targeting those GAM design areas are highly desirable for 
rationalizing further the causes underlying the discrepan¬ 
cies found between DFT-based and quantum chemistry 
methods [e.g., DMG and CCSD(T)], and in general for 
substantiating the role of many-body effects in the ad¬ 
sorption of gas molecules on atomically sparse environ¬ 
ments. In fact, many-body DFT methods themselves 
present also a number of remaining challenges like it can 
be the derivation of an analytic expression for the calcu¬ 
lation of atomic forces. (This ingredient is highly sought 
after for the realization of efficient geometry optimiza¬ 
tions and molecular dynamics simulations.) Work on this 
and other directions have been already initiated [49l |50| , 
thereby we expect that the application of many-body 
DFT methods will become routinary in the near future. 
We animate computational scientists to start consider¬ 
ing these advanced computational techniques in their 
planned studies of GAM. 


Concerning the generalization of DFT benchmark con¬ 
clusions attained in cluster systems to realistic materials, 
we have presented evidence for a number of related issues 
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like are the omission of long range dispersion interactions 
and t he d isguise of electronic band structure effects (see, 
Secs. IE and IV). An effective way of getting rid of those 
potential drawbacks is to recover the translational invari¬ 
ance in the simulations and to perform all the required 
calculations in a same periodic system. This possibility 
naturally points to the suite of exact ground-state quan¬ 
tum Monte Carlo methods (e.g., diffusion Monte Carlo 
-DMC-) as one of the most promising approaches for 
undertaking benchmark studies on the performance of 
DFT. In fact, the formalism of Bloch functions is already 
firmly established within QMC [SZllSH] and different im¬ 
plementations of this approach are available in several 
open-source packages (e.g., CASINO [59], QWALK [60] . 
and QMCPACK [M]). Also, the computational cost of 
QMC calculations is orders of magnitude more favorable 
than those of other popular quantum chemistry meth¬ 
ods. For these reasons, we envisage that QMC methods 
will play an increasingly more relevant role in prospective 
DFT benchmark studies of GAM. (Actually, the DMC 
method has already been applied to the study of metal 
hydrides [2281 1229L1287] and metal oxides |288H290j , two 
important families of materials within the context of hy¬ 
drogen storage and carbon capture.) It is worth notic¬ 
ing, however, that QMC methods are neither exempt of 
some important technical problems. These are essentially 
related to the tediousness found in the computation of 
atomic forces and the correction of numerical 

bias introduced by the nodal surface approximation [58] . 

An alternative to using periodic quantum chemistry 
methods, which is computationally more feasible but also 
tentative, may consist in designing model cluster sys¬ 
tems which mimic the electronic structure of targeted 
extended GAM (see Sec. IVBlI and where all the hi¬ 
erarchical calculations are subsequently performed. In 
following this approach, one should correct somehow for 
the electronic long-range exchange and correlations in the 
final adsorption energies. This could be done ad hoc by 
carrying out additional calculations in both periodic and 
cluster systems with a suitable first-principles method 
(e.g., FLRC functionals like a;B97X-D and BLYP-D3). 
Such a proposed long-range correction scheme is very 
much on the spirit of hybrid DFTiQCM approaches, 
which have been demonstrated to be successful on the 
simulation of extended systems (e.g., see Sec. IVB 2 and 
works [2851 1286] '). The reliability of this alternative 
benchmark approach, however, has not been yet fully 
assessed hence further work on this direction is certainly 
needed before claiming any progress. 


VI. GENERAL CONCLUSIONS 

In this critical review, we have presented abundant ev¬ 
idence showing that both electronic long-range exchange 
and long-range correlations play a decisive role in the 
adsorption of H 2 and CO 2 molecules on carbon-based 
materials and MOF. Even in situations where the cal¬ 


culated gas binding energies turn out to be large (i.e., 
Fbind >0.1 eV), the same conclusion holds to be true. 
Consequently, DFT exchange-correlation functionals em¬ 
ployed in the modeling of GAM must incorporate fea¬ 
tures which somehow correct for standard self-interaction 
errors and simultaneously reproduce dispersion forces. 
DFT Ej;c functionals of this class include the brand-new 
BLYP-D3, a;B97X-D, and vdW-optB88 variants, to cite 
a few. Other similar “full long range corrected” func¬ 
tionals can be tailored in principle by combining hybrid 
and dispersion correction schemes of one’s personal taste, 
producing so a suite of safe DFT approaches which can 
range widely on versatility and computational expense. 

As for standard DFT functionals (i.e., LDA and GCA), 
since these completely disregard electronic long-range ex¬ 
change and correlations, we recommend not to use them 
in the modeling of GAM and simulation of gas adsorption 
phenomena in general. In some particular circumstances 
LDA and GGA approaches may provide the correct quali¬ 
tative answers, however this is likely to occur fortuitously 
as a result of large error cancellations. Regrettably, cur¬ 
rent trends realized in modeling of GAM show that the 
use of LDA and GGA functionals is still quite widespread. 
We expect to motivate a change on this tendency with 
the present critical review. Also, we appeal to bring into 
new examination reported standard DFT predictions on 
GAM which are in conflict with observations. 

Concerning prospective DFT work, we animate re¬ 
searchers to consider the RPA-DFT and DFT-I-MBD 
methods in future studies of H 2 storage and carbon cap¬ 
ture materials. The reason for this is to substantiate with 
precision the role of many-body energy and Coulomb 
screening effects on the estimation of gas binding energies 
and equilibrium geometries. This is a completely new di¬ 
rection to take in modeling of GAM and, despite that 
the applicability of many-body DFT-based approaches 
remains yet limited, we believe that it will help in under¬ 
standing further the origins of the discrepancies found 
between DFT and quantum chemistry results. 

Finally, we recommend to be cautious in generalizing 
the conclusions of benchmark studies performed in model 
cluster systems to realistic GAM. In fact, the density of 
electronic states around the Fermi level and the HOMO 
and LUMO orbitals in chemically similar cluster and 
extended systems, even when considering medium size 
model complexes, can be very different. Those circum¬ 
stances easily can translate into completely different gas- 
GAM governing interactions. Also, electronic long-range 
interactions present in extended materials generally turn 
out to be disguised in model cluster systems. The safest 
strategy for avoiding these potential scaling problems is 
to consider quantum chemistry methods which, alike to 
DFT, can handle the simulation of periodic systems. In 
this context, quantum Monte Garlo emerge as one the 
most effective suite of exact ground-state methods for 
computation of gold-standard benchmarks in GAM. 
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